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Abstract. Results of optical and NIR spectral and photometric observations of a sample of candidate double- 
barred galaxies are presented. Velocity fields and velocity dispersion maps of stars and ionized gas, continuum 
and emission-line images were constructed from integral-field spectroscopy observations carried out at the 6m 
telescope (BTA) of SAO RAS, with the MPFS spectrograph and the scanning Fabry-Perot Interferometer. NGC 
2681 was also observed with the long-slit spectrograph of the BTA. Optical and NIR images were obtained at 
the BTA and at the 2.1m telescope (OAN, Mexico). High-resolution images were retrieved from the HST data 
archive. Morphological and kinematic features of all 13 sample objects are described in detail. Attention is focused 
on the interpretation of observed non-circular motions of gas and stars in circumnuclear (one kiloparsec-scale) 
regions. We have shown first of all that these motions are caused by the gravitational potential of a large-scale 
bar. NGC 3368 and NGC 3786 have nuclear bars only, their isophotal twist at larger radii being connected with 
the bright spiral arms. Three cases of inner polar disks in our sample (NGC 2681, NGC 3368 and NGC 5850) are 
considered. We found ionized-gas counter-rotation in the central kiloparsec of the lenticular galaxy NGC 3945. 
Seven galaxies (NGC 470, NGC 2273, NGC 2681, NGC 3945, NGC 5566, NGC 5905, and NGC 6951) have inner 
mini-disks nested in large-scale bars. Minispiral structures occur often in these nuclear disks. It is interesting that 
the majority of the observed, morphological and kinematical, features in the sample galaxies can be explained 
without the secondary bar hypothesis. Thus we suggest that a dynamically independent secondary bar is a rarer 
phenomenon than follows from isophotal analysis of the images only. 
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' 1. Introduction trarily oriented with respect to the primary (outer) one. In 
, the following years the isophotal analysis technique made 
l£j . More than a quarter of a century ago de Vaucouleurs it possib ie to detect secondary bars in optical (Wozniak et 
CS ; (1975) found that an inner bar-like structure was nested in aL; 19g5; Erwin & Sparke , 1999; Erwin & Sparke, 2003) 
the large-scale bar in the optical image of NGC 1291. The and NIR images of barred ga i axies (Friedli et ah, 1996; 
first systematic observational study of this phenomenon j ungwiert et aLj i997; Laine et ^ 2 002). Up to now, we 
was made by Buta & Crocker (1993), who published a list have found 71 candidate double-barred galaxies in the lit- 
of 13 galaxies, where the secondary (inner) bar was arbi- erature ( see references in Moiseev, 2001b). A new revised 
catalog of 67 barred galaxies including 50 double-barred 



Send offprint requests to: A.V. Moiseev e-mail: moisav@sao.ru ones is presented by Erwin (2004). Although various the- 

* Based on observations carried out at the 6m telescope of oretical studies exist, the secondary bar dynamics is still 

the Special Astrophysical Observatory of the Russian Academy f ar f rom being well understood. Shlosman et al. (1989) 

of Sciences, operated under the financial support of the Science haye shown that an additional circumnuclear bar may 

Department of Russia (registration number 01-43), at the be formed &g ft regult Qf instabilities of the gaseQUS ring 

2.1m telescope of the Observatorio Astrononico Nacional San T T . ., , , „ rr T r>\ • ^ 
_ . „,,:,»,. , r , , n ■ r , i m the Inner Lmdblad Resonance ILK within a large- 
Pedro Martir, Mexico, and from the data archive of the , , v ' ° 
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sometimes appear in stellar-gaseous simulations of galac- 
tic disks (Pfenninger & Norman, 1990; Friedli & Martinet, 
1993). New hydrodynamical simulations of the gas behav- 
ior in double bars, with an analysis of the problems arising 
during modeling and interpretation of the double-barred 
structures are presented by Maciejewski et al. (2002) and 
Shlosman & Heller (2002). 

Various observational data indicate that we are prob- 
ably seeing a new structural feature of barred galaxies in 
the case of double bars. However, the majority of obser- 
vational results has been based on photometric measure- 
ments only, when an elongated elliptical structure appears 
inside the large-scale bar. Moreover, using the isophotal 
analysis formalism (Wozniak et al., 1995) has allowed even 
"triple bars" (Friedli et al., 1996; Jungwiert et al., 1997; 
Erwin & Sparke, 1999) without any arguments about the 
dynamical behavior of such complex stellar systems. At 
the same time, the structures seen in direct images can 
be explained in a less exotic manner, without using the 
hypothesis of double or triple bars. Many structural fea- 
tures would be able to create an illusion of a secondary bar 
(Friedli et al., 1996; Moiseev, 2001b; Petitpas & Wilson, 
2002), such as for example, an oblate bulge, a complex 
distributions of dust and star formation regions in the 
circumnuclear region, or elliptical rings at the ILRs of a 
primary bar (Shaw et al., 1993). Therefore new observa- 
tional data are required to test our understanding of the 
secondary bar phenomenon. 

Kinematic data (the stellar and gas line-of-sight veloc- 
ities and velocity dispersion) over the circumnuclear re- 
gions of barred galaxies are of crucial importance for this 
kind of studies. Recently, Emsellem et al. (2001a) have 
studied the kinematics of double-barred galaxy candidates 
and reach a conclusion on dynamical decoupling of the 
circumnuclear regions in NGC 1097, NGC 1808 and NGC 
5728. This conclusion is based on the fact that peaks of 
the relative line-of-sight velocities and stellar velocity dis- 
persion drops are seen in the circumnuclear regions from 
long-slit cross-sections. Wozniak (2000) has supposed that 
"counter-rotation" of stars, observed in the nuclear region 
of NGC 5728, may be associated with secondary bar in- 
fluence on stellar kinematics. All these features observed 
in the onc-dimcnsional data may be explained without 
using the hypothesis of secondary bars. This may be re- 
lated to the specific mass distributions in the inner kilo- 
parsec and to non-circular motions of stars and gas within 
the primary bar only (a^-families of stellar orbits lead to 
"counter-rotation" properties). Also Wozniak et al. (2003) 
in their self-consistent N-body simulations have shown 
that the velocity dispersion drops observed by Emsellem 
et al. (2001a) can be independent on the presence of a sec- 
ondary bar and are reproduced in the model of a single- 
barred galaxy. 

Since the motions of stars and gas inside the bar re- 
gion are strongly non-circular, and such objects are non- 
axisymmetric by definition, panoramic (also named 2D, 
3D or integral-field) spectroscopy should give us more de- 
tailed information on circumnuclear kinematics in com- 



pared to "classical" (long-slit) spectroscopy. Integral-field 
spectroscopy makes it possible to map two-dimensional 
fields of line-of-sight velocities and velocity dispersion. 
Therefore, during 2000-2002 we have carried out a series 
of observational programs on 2D-kinematics and morphol- 
ogy of double-barred galaxies candidates. The main goal 
was to find an answer to the question - Are the secondary 
bars dynamically decoupled systems? 

This problem is complicated and ambiguous because of 
the lack of theoretical studies, where the "observed" ve- 
locity fields in the double-barred galaxies had been mod- 
eled. It is important to collect a full sample of candidate 
double-barred galaxies with observable 2D-kinematics of 
stars and gas. These observations must be used to search 
for common dynamical (kinematic) features in these kinds 
of galaxies. Recently, Shlosman, & Heller (2002) have writ- 
ten: "..clearly, the most promising method in detecting the 
nuclear bars is two-dimensional spectroscopy of the central 
kiloparsec, which can reveal the underlying kinematics" . 

Our observational program was finished in 2002. 
Independently similar observations were started at the 
3.6m CFHT using the integral-field spectrographs OASIS 
and TIGER. Recently, the velocity fields of stars and ion- 
ized gas in some double-barred candidates have been pre- 
sented by Emsellem & Friedli (2000), Emsellem et al. 
(2001b) and Emsellem (2002). The velocity fields of the 
molecular gas were also investigated in double-barred can- 
didates NGC 2273 and NGC 5728 (Petitpas & Wilson, 
2002) and NGC 4303 (Schinnerer et al., 2002). 

In this work, we present the observational data 
(integral- field spectroscopy and NIR/optical surface pho- 
tometry) for 13 galaxies and discuss the problem of central 
region dynamical decoupling. The paper is structured as 
follows. In Sect. 2 we describe the observations and data 
reduction technique. The methods of analysis of the veloc- 
ity fields and surface brightness distributions are consid- 
ered in Sect. 3. In Sect. 4 we present the notes for individ- 
ual observed galaxies and discuss their common properties 
in Sect. 5; a short conclusion is drawn in Sect. 6. 

2. Observations and data reduction 

We have compiled our sample from Moiseev (2001b) by 
applying the following restrictions: S > 0, and the size of 
the major axis of the probable secondary bar is smaller 
than or equal to the MPFS field of view (16" x 15"). We 
have obtained full observational data for 13 galaxies, that 
is for about one third of the candidates in the Northern 
sky. 

2.1. Integral-field spectroscopy with the MPFS 

Circumnuclear regions of all the galaxies were ob- 
served at the 6m telescope of the Special Astrophysical 
Observatory, Russian Academy of Sciences (SAO RAS) 
with the Multipupil Integral Fiber Spectrograph (MPFS, 
Afanasiev et al., 2001). A description of the spectro- 
graph is also available through the Internet at SAO WEB- 
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page http : //www. sao . ru/hq/lsf vo/devices .html. The 

spectrograph takes simultaneous spectra from 240 spatial 
elements (constructed in the form of square lenses) that 
form an array of 16 x 15 elements on the sky. The field of 
view is 16" x 15", and the angular size is l"/element. 

Three galaxies, NGC 2273, NGC 3368, and NGC 3945 
were observed in two different positions of the multilcns 
array on the sky, with each frame containing a galactic nu- 
cleus. After preliminary data reduction, the "data cubes" 
were combined. The resulting field of view was 16" x 18" 
for NGC 2273, 18" x 18" for NGC 3368 and 23" x 23" for 
NGC 3945. Together with target spectra, we took a night- 
sky spectrum from an area located 4'5 from the center of 
the field of view. The detector was a TK1024 1024 x 1024 
pixels CCD array. The spectrograph reciprocal dispersion 
was 1.35A/pixel and the spectral resolution was about 4A. 
The log of MPFS observations is given in Table 1 . It con- 
tains the names of galaxies, dates of the observations, the 
total exposures T exp , the spectral range A A and the seeing 
quality during the exposures. 

We reduced the observations by using the software de- 
veloped at the SAO RAS by V.L. Afanasiev and running 
in the IDL environment. The primary reduction included 
bias subtraction, flat-fielding, cosmic-ray hits removal, ex- 
traction of individual spectra from the CCD frames, and 
their wavelength calibration using a spectrum of a He- 
Ne-Ar lamp. Subsequently, we subtracted the night-sky 
spectrum from the galaxy spectra. The spectra of spec- 
trophotometric standard stars were used to convert the 
fluxes into absolute energies. 

The reduced spectra were represented as "data cubes" : 
each spatial element of the two-dimensional field has an in- 
dividual 1024-channel spectrum. We constructed the maps 
of line intensity and light-of-sight velocity fields in the H3 , 
[O III] A4959, 5007 A and/or [N II]A6548, 6583 A lines by 
means of the Gaussian fitting of the emission line profiles. 
A double-gaussian model was used for doublet separa- 
tion. The absolute accuracy of the velocity determination, 
evaluated from the air-glow lines wavelengths, was about 
10 — 15kms _1 . For the continuum map construction we 
summed the fluxes in the spectral ranges free of emission 
lines (5600 - 5800 A for "green" and 6350 - 6450A for 
"red" spectra). 

The line-of-sight velocity and velocity dispersion fields 
for the stellar component were constructed by means of a 
cross-correlation technique, modified for the MPFS data 
(Moiseev, 2001a). The used spectral region included some 
stellar absorption features; MgIA5175 A, FeIA5229 A, 
FcI+CaIA5270 A, and NaIA5893 A among others. The 
spectra of G8-K3 giant stars and the twilight sky, ob- 
served during the same nights, were obtained as velocity 
templates. To take into account the variations of the in- 
strumental contour (see Moiseev, 2001a)the observations 
of the template, the images of stars were out-of-focus 
to fill the MPFS field of view. The accuracies were of 
~ lOkrns -1 , and 10 — 20kms~ 1 respectively for the ve- 
locity and velocity dispersion determinations (the errors 
are dependent on the cross-correlation peak amplitudes). 



2.2. Observations with the Fabry-Perot Interferometer 

Six galaxies with bright emission lines were observed 
at the 6m telescope with the scanning Fabry-Perot 
Interferometer (IFP). The Queensgate interferometer ET- 
50 was used in 235th interference order (for the H Q line). 
The free spectral range between neighboring orders (in- 
terfringe) was about 28 A (1270 krns" 1 ). For preliminary 
monochromatization, a set of narrow-band filters, with 
a FWHM = 12 — 18 A and centered on the spectral 
region containing rcdshifted H a or [N IIJA6583 A emis- 
sion lines were used. During the observations we succes- 
sively take 32 interferometric images of an object with 
different gaps between IFP plates. The spectral channels 
have a width of SX ~ 0.9A (~ 40kms _1 ), the spec- 
tral resolution (the width of instrumental contours) was 
FWHM w 2.5 A (~ 110 krns- 1 ). The detector was a 
CCD TK1024 (1024 x 1024 pixels). The CCD was oper- 
ated with 2x2 binning for reading-out time economy. 
Therefore each spectral channel has a 512 x 512 pixel for- 
mat. 

In February and March, 2000, the IFP was installed 
into the parallel beam, inside the focal reducer (Dodonov 
et al., 1995; Moiseev, 2000). The resulting focal ratio at 
the prime focus of 6m telescope was F/2A, with a field of 
view of 5'8 and a scale of 0.68"/px. 

In September and November 2000, the ob- 
servations were carried out with the new multi- 
mode focal reducer SCORPIO. A brief de- 
scription of this device is given on Internet 
(http : //www . sao . ru/hq/moisav/scorpio/ scorpio . html), 
also the IFP mode in SCORPIO is described by Moiseev 
(2002a). With this configuration, the focal ratio was 
F/2.9, with a field of view of 4.'8 and a scale of 0.56"/px. 
The log of the IFP observations is presented in Table 2. 

To reduce the interferometric observations we used a 
custom development software (Moiseev, 2002a), running 
in the IDL environment. After primary reduction (bias, 
flat-fielding etc.), removing the night-sky emission lines 
and converting to the wavelength scale, the data were pre- 
sented as "data cubes" . In such a cube, to every pixel of 
the 512 x 512 field, a 32-channel spectrum is attached. The 
"data cubes" were smoothed by Gaussians with a width, 
in the spectral domain, of FWHM — 1.5 channel, and 
in the spatial (sky-plane) domain of FWHM = 2 — 3 
pixels under the ADHOC package 1 . The velocity fields of 
the ionized gas, and the images in the emission line were 
mapped by means of a Gaussian fitting of the emission line 
profiles. Moreover, we created the images of the galaxies 
in the "red" continuum near to the emission lines. 

2.3. NIR photometry at 2.1m telescope 

The NIR surface photometry of the sample galaxies was 
carried out with the CAMILA camera (Cruz-Gonzalez et 

1 ADHOC software was written by J. 
Boulesteix (Observatoire de Marseille). See 
http : //www-obs . enrs-mrs . f r/ ADHOC/ adhoc .html) 
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Table 1. Log of MPFS observations. 
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13 /05 /2002 


3 x 1200 


5900 
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1.5' 


NGC 4736 


11/05/2000 


3 x 1200 


4840 


-6210 


1.8' 


NGC 5566 


10/05/2000 


3 x 1200 


4840 


-6210 


2.3' 


NGC 5850 


13/08/2001 


3 x 1200 


5710 


- 7080 


2.7' 




13/05/2002 


3 x 1200 


4700 


- 6070 


1.2' 


NGC 5905 


30/03/2000 


8 x 1200 


4840 


-6210 


2.0' 


NGC 6951 


03/09/2000 


5 x 1200 


4840 


-6210 


2.2' 


NGC 7743 


02/09/2000 


3 x 1200 


4840 


-6210 


2.0' 



Table 3. Log of CAMILA observations. 



Name 


Date 


Mode 


Filters 


Seeing 


NGC 2273 


09/03/2001 


(F/4.5) 


J, H, K' 


3.5" 


NGC 2681 


09/03/2001 


(F/4.5) 


J, H, K' 


2.8" 




18/03/2000 


(F/13.5) 


J, H, K' 


1.5" 


NGC 2950 


12/03/2001 


(F/4.5) 


J, H 


3.0" 




18/03/2000 


(F/13.5) 


J, H, K' 


1.5" 


NGC 3368 


12/03/2001 


(F/4.5) 


J, H 


2.2" 




18/03/2000 


(F/13.5) 


J, H, K' 


1.5" 


NGC 3786 


19/03/2000 


(F/13.5) 


J, H 


1.5" 


NGC 3945 


09/03/2001 


(F/4.5) 


J, H 


4.0" 




18/03/2000 


(F/13.5) 


J, H, K' 


1.5" 


NGC 4736 


20/03/2000 


(F/4.5) 


J, H, K' 


1.5" 


NGC 5566 


18/03/2000 


(F/13.5) 


J, H, K' 


1.5" 




12/03/2001 


(F/4.5) 


J, H 


1.6" 


NGC 5850 


18/03/2000 


(F/13.5) 


J, H, K' 


1.5" 


NGC 5905 


19/03/2000 


(F/13.5) 


J, H 


2.5" 



Table 2. Log of IFP observations. 



Name Date 


Texp-, SGC 


Line 


Seeing 


NGC 470 02/11/2000 


32 x 200 


H a 


2.1" 


NGC 2273 03/11/2000 


32 x 200 


H a 


2.0" 


NGC 3368 28/02/2000 


32 x 150 


H a 


2.7" 


28/02/2000 


32 x 200 


[N II] 


3.5" 


NGC 3945 03/11/2000 


32 x 190 


[Nil] 


2.0" 


NGC 4736 02/03/2000 


32 x 150 


H a 


2.7" 


02/03/2000 


32 x 180 


[N II] 


3.5" 


NGC 6951 23/09/2000 


32 x 120 


H a 


2.5 


23/09/2000 


32 x 120 


[N II] 


2.5" 



al., 1994) in the JHK' bands, on the 2.1m telescope of the 
Observatorio Astronomico Nacional, San Pedro Martir, 
Mexico. CAMILA is equipped with a NICMOS3 256 x 256 
pixels detector. Observations were carried out with two 
telescope focuses, with different focal ratios. In the F/13.5 
mode the field of view was 1.'3 and the pixel size 0730. The 
F/4.5 mode provided a field of view of 3.'6 with a scale 
of O'.'85/px. The log of the NIR observations is given in 
Table 3. 

Each observation consisted of a sequence of object and 
sky exposures, with the integration time of an individual 
exposure limited by the background level, which was kept 
well below the non-linear regime of the control electronics. 
The nearby sky frames were taken with the same exposure 
times and at adjacent positions to the object, through 
an observing sequence that alternates the object and sky 
frames. 

The routines from Image Reduction and Analysis 
Facility (IRAF) were used in the reduction and analysis 
of all the data. The image processing involved subtraction 
of sky frames, division by flat field frames, registration of 
the images to a common coordinate system and stacking 
all the images in a filter. Bias subtraction was carried out 



during the data acquisition. Some of the images had a hor- 
izontal band at the joints of the individual quadrants of 
the NICMOS 3 array detector. We eliminated this feature 
by subtracting a 255x1 pixels median smoothed image 
from the original images. 

2.4. Direct imaging on the 6m telescope 

Optical direct imaging observations of some galaxies were 
carried out with the focal reducer SCORPIO (Sect. 2.2) at 
the 6m telescope. The detector (CCD TK1024, 1024 x 1024 
pixels) provided a field of view of 4.'8 with a pixel scale 
of 0'.'28. The log of the observations is presented in 
Table 4. We used 3 filters: V (Johnson), R c (Cousins) and a 
medium-band red filter (SED755) for the continuum imag- 
ing of the objects (without any pollution from the emis- 
sion lines of the object). The SED755 filter has a width 
of FWHM = 220 A and is centered at 7560 A. The data 
reduction included bias and flat-field corrections, and cos- 
mic hits removing. 

2.5. Hubble Space Telescope data 

We also used the high-resolution images of the sample 
galaxies available from the HST Archive. The HST images 
were obtained with the optical WFPC2 and the infrared 
NICMOS cameras. WFPC2 has an image scale of 0.1"/px 
and a field of view of 2'.7 and 0.045/px for the Planetary 
Camera (PC), with a field of view of 36". The NICMOS 
detector includes three chips. Its pixel scale and field of 
view are 0//043 and 11" for NIC1, 07075 and 19" for NIC2, 
and 072 and 51" for NIC3. The information about the 
archival images is given in Table 5. Here ID is an archive 
identification of the observational program. 
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Table 4. Log of direct imaging observations at the 6m tele- 
scope 



Name Date Filters Seeing 



NGC 470 03/11/2000 


V, 


R c , SED755 


1.4" 


NGC 2273 03/11/2000 


V, 


Re 


1.0" 


NGC 2681 25/10/2000 


V 




1.0" 


NGC 2950 04/11/2000 


v, 


Re 


1.5" 


NGC 3945 03/11/2000 


v, 


R c 


1.0" 


NGC 6951 04/11/2000 


v, 


R c , SED755 


1.2" 


NGC 7743 03/11/2000 


v, 


Rc , SED755 


1.5" 



Table 5. HST high resolution direct images 



Name 


Date 


Camera 


Filters 


ID 


NGC 


2273 


05/02/1997 


WFPC2 


F791W 


6419 






03/04/1998 


NIC2 


1 160W 


7172 






03/04/1998 


NIC3 


F164N, F166N 


7172 


NGC 


2681 


07/06/1998 


NIC3 


F187N, F160W 


7919 






03/09/2000 


WFPC2 


F300W 


8632 


NGC 


2950 


26/01/1999 


WFPC2 


F450W, F555W, 


6633 










F814W 




NGC 


3368 


04/05/1998 


NIC2 


F160W 


7331 






08/05/1998 


NIC2 


F110W 


7331 






12/12/2000 


WFPC2 


F814W 


8602 


NGC 


3786 


30/03/1995 


WFPC2 


F606W 


5479 






29/04/1998 


NIC1 


F110W, F160W 


7867 


NGC 


3945 


11/05/1997 


WFPC2 


F450W, F555W 


6633 










F814W 




NGC 


4736 


02/07/2001 


WFPC2 


F814W 


9042 


NGC 


5566 


04/04/1999 


WFPC2 


F606W 


6359 


NGC 


6951 


30/10/2000 


WFPC2 


F814W 


8602 






18/11/2000 


WFPC2 


F606W 


8597 


NGC 


7743 


30/05/1994 


WFPC2 


F547M 


5419 






11/10/1994 


WFPC2 


F606W 


5479 






16/09/1997 


NIC2 


F160W 


7330 



2.6. Long-slit observations of NGC 2681 

When all the data were reduced and analysed we con- 
cluded that the photometric and MPFS data were not 
enough for understanding the outer morphology in NGC 
2681; first of all, the PA of the disk was very uncertain 
(Sect. 4.4). For this goal additional observations of this 
galaxy were performed in April-May of 2003 at the 6m 
telescope with the SCORPIO focal reducer in a long-slit 
mode. The detector was the EEV42-40 2048 x 2048 pixels 
CCD array. The spectrograph reciprocal dispersion was 
0.85A/pixel, and the spectral resolution was about 6A. 
The slit of 1.2" x 5.3' was centered onto the nucleus under 
three fixed orientations near possible minor (PA=30° and 
61°) and major (PA=121°) axes. The log of the observa- 
tions is given in Table 6. Though the seeing was 2.5" — 3.5" 
on April 7, this is not important for the study of the kine- 
matics of the outer regions.. The primary reduction in- 
cluded bias subtraction, flat-fielding and cosmic-ray hits 
removal. Then the individual spectra were sampled along 
the slit with a bin of 0.6" at distances r < 10" from the 



Table 6. Log of long- slit observations 



Name 


Date 


PA 


T eX pi SGC 


Seeing 


NGC 2681 


07/04/2003 


30° 


4 x 600 


2.6" 


NGC 2681 


07/04/2003 


61° 


3 x 900 


3.5" 


NGC 2681 


01/05/2003 


121° 


3 x 1200 


1.5" 



center and with a bin of 2.5" at larger distances. The spec- 
trum of a He-Ne-Ar lamp was used for wavelength cali- 
bration. The line-of-sight velocity distributions along the 
slit were constructed by means of a cross-correlation tech- 
nique. The spectra of the K0III stars observed in the same 
nights were used as velocity templates. The errors of the 
velocity measurements are from lOkms -1 in the galactic 
center to 20 — 30kms _1 at distances of r = 40 — 45". 

3. Analysis of the velocity fields and photometric 
data 

The isophotc analysis of the images was performed with 
the FITELL program written by V.V. Vlasyuk (SAO 
RAS) which uses the well-known algorithm of Bender & 
Moellenhoff (1987). It yields a surface brightness value (/) 
and the orientation parameters of elliptical isophotes: po- 
sition angle PA, and ellipticity e = 1 — b/a, where a and b 
are the semi-major and semi-minor axis of ellipses respec- 
tively. 

The velocity fields were analyzed by means of the 
"tilted-ring" method (Bcgcman, 1989), with fixed sys- 
temic velocity and position of the center as described 
by Moiseev & Mustsevoi (2000). Here we briefly describe 
the method. The velocity field is broken up into elliptical 
rings, with a 1 — 1.5" width aligned with the direction of 
the line of nodes. Using the x 2 minimization of the devia- 
tions of the observational points from the model, we calcu- 
late all 6 parameters characterizing pure circular rotation 
at a radius R in the galactic plane. They are: the coor- 
dinates of the rotation center on the sky plane (xq, yo), 
the systemic velocity V sys , the position angle of the major 
axis PAdyn (dynamical position angle), the galactic plane 
inclination to the linc-of-sight i, and the mean rotation ve- 
locity V ro t- However, since all galaxies contain large-scale 
bars in the inner and spiral arms in the outer regions, the 
"circular" velocity fields in the observed galaxies must be 
distorted by non-circular motions. The non-circular mo- 
tions cause in the systematic errors when determining the 
disk orientation parameters (see Appendix in Lyakhovich 
et al., 1997; Barnes & Sellwood, 2003; Fridman et al., 
2004). 

In such a situation it is convenient to fix some 
model parameters. A good first approach is to assume 
xq, yo, V sys , andi = const, leaving only two free param- 
eters for fitting in each ring: V rot and PAdyn ■ The 
PAdyn deviation from the line of nodes of the whole disk 
(PAq) is used to characterize the type of non-circular mo- 
tions (oval distortion, polar disk, etc.). In the case of the 
motions of gas clouds in a bar potential, the observable 
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PAdyn ceases to be aligned with the disk's line of nodes 
(Chevalier & Furenlid, 1978, Moiseev & Mustsevoi, 2000). 
A turn of the "dynamical axis" must occur in the oppo- 
site direction, with respect to position angle of the in- 
ner isophotes, or to the "photometric axis" (Moiseev & 
Mustsevoi, 2000). The same effects should be also observed 
in the velocity field of the stars, as is shown by numerical 
simulations (Miller & Smith, 1979; Vauterin & Dejonghe, 
1997), and by analytical calculations of the stellar mo- 
tions under the effects of a triaxial gravitational potential 
(Monnet et al., 1992). On the other hand, in the case of a 
pure circular rotation in a warped (polar, inclined) inner 
disk, the dynamical axis must follow the photometric one. 

In determining the set of optimal galactic disk orien- 
tation parameters we were faced with some problems. For 
the systemic velocity V sys we took its mean value over all 
radii. As a first approach, for the center position x ,yo 
we took the coordinates of the symmetry center of the 
velocity field. If this "dynamical center" is at a distance 
less than 1 — 1.5" from the center of the inner isophotes 
of the continuum image, then the photometric center is 
accepted as the ultimate position of the center. In NGC 
470 (Sect. 4.2), the difference between the positions of the 
dynamical and photometric centers exceeded the error of 
their measurement (the case of a lopsided galaxy). For this 
galaxy we adopted the dynamical center as the center of 
rotation. 

The main problem is to accurately calculate the val- 
ues of the position angle of the line of nodes {PA ) and 
the inclination (i). To use the averaged parameters of the 
outer isophotes is not always correct for this, because the 
outer parts of the spiral structure distort these isophotes in 
galaxies possessing well-developed spiral arms. The choice 
of radius range for which the average values of PA and e 
are calculated is also not obvious. In addition, some galax- 
ies (NGC 3368, NGC 3945) have ring-like structures at 
large distances from the center. If a ring has a resonant ori- 
gin, then it could be elliptical in the galactic plane (Buta, 
1996) with uncertain ellipticity for an observer. Finally, 
the errors in flat-fielding and background-removing proce- 
dures (mainly for JHK' images) also distort the shape of 
the low-brightness isophotes of the outer regions of galax- 
ies. Therefore, we used the "spiral criterion", proposed by 
Fridman et al. (2004), for PA and i determinations for 
the main part of our sample. This method is based on ap- 
plying of a Fourier-series expansion to the surface bright- 
ness distribution in the galactic plane I(R, <p), where R is 
the radial distance, and ip the azimuthal angle: 

N 

I{R,<p) = A {R)+Y^ A m (R) cos(m (p + 4> m (R)) (1) 

m— 1 

A m , and (f> m are the amplitude and phase of the m-th har- 
monic, and N — 8 to 12 is the maximum number of the 
terms. If a m-armed spiral structure appears on the image, 
then the line of the maximum of the harmonic number m 
should also be spiral. Clearly, in the bar region the line 
of the maximum of the m = 2 harmonic must be straight 



Table 7. Parameters of disk orientation 



Name 


PA 




i, 


n 


NGC 470 


157 


±i 


55 


±i 


NGC 


2273 


58 


±2 


50 


±2 


NGC 


2681 


148 


±5 


25 


±5 


NGC 


2950 


119 


±1 


50 


±2 


NGC 


3368 


135 


±5 


48 


±3 


NGC 


3786 


247 


±5 


61 


±2 


NGC 


3945 


157 


±2 


52 


±2 


NGC 


4736 


298 


±6 


31 1 




NGC 


5566 


210 


±3 


64 


±2 


NGC 


5850 


335 


±10 


37 


±5 


NGC 


5905 


130 2 




40 2 




NGC 


6951 


140 


±6 


41 


±4 


NGC 


7743 


270 


±5 


40 


±3 



1 Mulder (1995) 

2 van Moorsel (1982) 



and elongated along the position angle of the bar. It ap- 
pears that the shape of the spiral, traced by the maximum 
of the harmonics, is very sensitive to the selection of the 
orientation parameters. It can be shown that, if for a two- 
armed galaxy the inclination i is underestimated, then the 
line of the maximum of the rn = 2 harmonic at large radii 
will be fixed near the direction of the line of nodes. 

Similarly, if the adopted PAq differs from the "true" 
value, then the line of the maximum of the m — 2 har- 
monic will be approximate to the direction of the "true" 
PAq in the outer part of the disk. Our experience shows 
that for moderate inclinations (i = 30 — 70°), the pa- 
rameters of disk orientations may be calculated with an 
uncertainty of a few degrees. A more detailed description 
of this method is given in Fridman et al. (2004). In Table 7 
we present the PA and i for our sample galaxies. During 
the search for the orientation parameters we tried to ob- 
tain some consistency between all three methods (isopho- 
tal analysis, analysis of the IFP velocity fields, and the 
"spiral criterion"). Nevertheless, preference was given to 
the "spiral criterion" . The details of the orientation results 
are considered in the discussion of individual galaxies. 

4. The results 

4.1. The Atlas 

The images and the results of the velocity field analysis 
arc presented in Figs. 1-18. The left panels of Figs. 1-13 
show the MPFS data. The top row contains the contin- 
uum images, the velocity field of the stars, and the map 
of the line-of-sight velocity dispersion (<r*). The image in 
the brightest emission line of the ionized gas, the velocity 
field of this emission line and the map of the velocity dis- 
persion of the ionized gas {cr gas ) are given in the middle 
row. For the last map the width of the instrumental con- 
tour was taken into account. There are no ionized maps 
for NGC 2950, NGC 3368, NGC 4736, and NGC 5566, 
because emission lines are absent in the MPFS spectra of 
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these galaxies. A gray scale is drawn, in kms -1 , for the 
velocity fields and for the velocity dispersion maps; while 
it is in ergss -1 cm _2 /D" for the continuum and emission 
line images, except for NGC 5850 in the [N II] emission 
lines, where the scale is in arbitrary counts. The cross 
marks the position of the adopted center of rotation in all 
MPFS maps. At the bottom we show the radial depen- 
dences of the rotation velocity V rot and the position angle 
on the dynamical axis PAd yn for the stellar and gaseous 
components. The dotted line marks the adopted position 
angle of the line of nodes (PA ). 

Examples of the images and isophotal analysis results 
are shown in the right panels of Figs. 1-13. The left col- 
umn contains the ground-based images with the largest 
field of view, while the right column contains the same 
data for the circumnuclear regions (r < 5 — 20"), observed 
with the HST. The best-resolution ground-based images 
are presented for those galaxies where there are no HST 
data. The logarithm-scaled images with the correspond- 
ing isophotes are presented in the top row of the panels. 
In the middle row we present the residual brightness distri- 
butions, after the removal of elliptical isophotes from the 
top images. In this case, the gray scale is linear. The thick 
dotted line shows the maximum of the m = 2 harmonic 
for the surface brightness, overlapped onto the ground- 
based images. At the bottom, we plot the parameters of 
the isophotes: PA and ellipticity versus major semi-axis. 
The scatter of the measurements, through the different 
filters, gives us some information about the errors of the 
elliptical model. The horizontal dotted line corresponds 
to the PAq value, and the vertical dotted line (in the left 
plot) marks the radial range, which is zoomed in on the 
right-hand plots. 

The results of the IFP observations are shown in 
Figs. 14-18, which include the continuum and H Q or/and 
[N II] emission line images. A square-root gray-scale is 
used. The cross marks the dynamical center. Also shown 
is the velocity field with isovelocities. The thick black con- 
tour corresponds to the systemic velocity, while the thin 
contours indicate line-of-sight velocities in steps of ±50 
kms -1 . The last plot shows the rotation curve (V ro t) and 
the kinematic position angles (PAd yn ) ■ Here the dotted 
line marks the PA value. For NGC 3368, NGC 4736, and 
NGC 6951 we present data in the H Q and [N II] emission 
lines. 

It is necessary to note briefly the errors in the V ro t de- 
termination. The formal errors of the model usually arc 
1 — 5kms _1 because every elliptical ring for fitting con- 
tains from several tens up to hundreds of points (the latter 
for IFP velocity fields). Nevertheless, the rotation curve, 
determined in the frame of pure circular motions, will be 
affected by systematic errors on the order of the spread 
in the observational velocities around the average values 
at the given radius (Lyakhovich et al., 1997; Fridman ct 
al., 2004). The error bars of V rot (r) in our figures must be 
interpreted as the deviation of the line-of-sight velocities 
from the models of circular rotation in each ring; i.e. the 
mean value of non-circular velocities. Therefore the "real" 



circular rotation curve should coincide with our curves 
within the errors which are given in Figs. 14-18. 

In the following subsections we will present the results 
for individual sample galaxies. 

4.2. NGC 470 

The adopted orientation parameters, PA = 149 ± 3°, 
i = 51 ± 4°, have negligible differences from the measure- 
ments of Garcia-Gomez & Athanassoula (1991). However, 
the parameters from Table 7 allow us to obtain a more 
smooth well-ordered spiral of the maximum of the m = 2 
harmonic. In the circumnuclear region, PA of the SED755 
isophotes matches the NIR data of Friedli et al. (1996), but 
it deviates by more than 25° from the measurements in the 
V and R c bands. This can be related to the dust extinc- 
tion, which is more significant in the optical bands. The 
deviations of PAd yn from the line of nodes at r = 8 — 20" 
(from the IFP data, see Fig. 14) arc explained by us as an 
effect of the large-scale (primary) bar on the ionized gas 
velocities. The turn of PAd yn occurs in the opposite direc- 
tion with respect to the inner isophotes PA (Fig. 1), in 
agreement with existing models of the gas flowing within 
the bars (see Sect. 3). Also the MPFS data reveal that 
the dynamical axis, for the stellar component, turns in 
the same direction as the gaseous (in H a ) one at r > 5". 

The center of rotation (dynamical center) of the MPFS 
velocity fields of stars and gas is shifted by ~ 4" (~ 0.6 
kpc), in the NW direction, from the photometric center. 
The discrepancy between the positions of both centers is 
considerably larger than the measurement errors or a pos- 
sible effect of dust absorption. The difference in the lo- 
cations of the centers is also observed in the large-scale 
H a velocity field, which confirms the result of Emscllem 
(2002) . Such a discrepancy in the photometric and dynam- 
ical centers locations or between the centers of the inner 
and outer isophotes (lopsided-structure) have been seen in 
numerous galaxies (Richter & Sancisi, 1994). The nature 
of the lopsidedness could be different at different radial 
scales, from the circumnuclear disk precessing around a 
super-massive black hole to development of a spiral az- 
imuthal disturbance m = 1 mode in the large-scale stellar- 
gaseous disk (Junqucra & Combes, 1996; Emscllem et al., 
2001a; Emsellem, 2002). 

Turnbull et al. (1997) have proved that NGC 474 shows 
signs of interaction with NGC 470. They also noted that 
both galaxies are perhaps embedded in a common HI- 
cloud. Hence, the excitation of the non-axisymmetrical 
harmonics in NGC 470 may be caused by a tidal inter- 
action with the companion. 

An alternative explanation can also be proposed for 
the lopsidedness since the distorted region is relative small 
and located inside the large-scale bar. It is possible to sup- 
pose an abnormal development of the non-axisymmetrical 
harmonics in the surface brightness during the secular evo- 
lution of a barred galaxy. A similar shift between the 
isophotal center and the dynamical center is observed 
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Fig. 1. (To be seen in landscape)NGC 470. The left panel shows the results of the MPFS-observations. Here the top row 
contains the continuum images, the velocity field of the stars, and the map of their velocity dispersion. In the middle row the 
emission line image, the velocity field of the ionized gas and the map of the velocity dispersion of the ionized gas are given. At 
the bottom the radial dependences of the rotation velocity and the PA of the dynamical axis are shown. The dotted line marks 
the line of nodes (PAo). In the right panel the left column contains the large-scale image, while the right column contains the 
image of the nuclear region. The logarithm-scaled images with the corresponding isophotes are presented in the top row of the 
panel. The residual brightness distributions are plotted in the middle row. The thick dotted line shows the maximum of the 
m=2 harmonic. At the bottom the PA and ellipticity of the isophotes are plotted. The horizontal dotted line corresponds to 
the PAo value, and the vertical dotted line (in the left plot) marks the radial range, which is zoomed in on the right-hand plots. 

Fig. 2. The same as Fig. 1 for NGC 2273 (To be seen in landscape) 

Fig. 3. The same as Fig. 1 for NGC 2681 (To be seen in landscape) 



in the central regions of some barred galaxies (Zasov & 
Khopcrskov, 2002). 

Wozniak et al. (1995) assumed that the isophotal turn 
at r < 7" was connected with a secondary bar. However 
later Friedli et al. (1996) have explained this feature by 
a triaxial bulge in the primary bar. This explanation is 
problematic from our point of view, because the position 
angle of the inner isophotes (at the r = 2 — 6") is close 
to the disk PA value. Also the amplitude of isophotal 
PA deviations from PAq decreases sharply when pass- 
ing from short wavelengths (V) to longer ones (SED755). 
Such behaviour of the isophotes can be explained by dust 
extinction. Moreover, the apparent ellipticity of the inner 
isophotes e — 0.3 — 0.45 (Fig. 1) corresponds to the projec- 
tion of a circular disk onto the sky plane (if problems with 
bulge and outer bar deprojection are taken into account). 
In Sect. 5.3 we decompose the surface brightness distri- 
bution in NGC 470 and justify the assumption. Thus we 
think that there is a circumnuclcar disk within the large- 
scale bar (see Sect. 5.3 for additional arguments). In the 
H Q velocity field the turn of PAd yn occurs in the same di- 
rection with respect to the PA of isophotes at r < 5 — 7" 
(Fig. 14). Such a dynamical feature is also in agreement 
with circular motions in the decoupled disk which may be 
slightly inclined to the galactic plane (Sect. 3). An asym- 
metrical m = 1 harmonic is developed in this disk, which 
also has a large concentration of molecular CO (Sofue et 
al., 1993), as well as of ionized hydrogen (see the H a - 
image on Fig. 14). The decoupling of the central disk may 
be caused by its location at the ILR of the large-scale bar. 

Vega Beltrah et al., (2001) gave an estimation of the 
central velocity dispersion of stars; er* = 56 ± 31kms _1 
from their long-slit observations of NGC 470. However, 
such a low value of cr* is doubtful, as it is difficult to 
explain the fact that the rotation velocity (from our 
MPFS data) of stars in the central region is smaller by 
~ 50 kms -1 than the gas velocity. Our measurements sup- 
port a central velocity dispersion of er* w 140—150 km s~ 1 , 
which looks more realistic and agrees with estimations of 
Heraiideau et al. (2001). The disagreement with the Vega 
Beltrah et al. (2001) results may be due to a template mis- 
matching problem, at was already marked in their article. 



In the H Q velocity field, at distances of r = 12 — 20" to 
the East from the center, we observe a region with line- 
of-sight non-circular motions about 50 — 60kms~ 1 . The 
profile of the H Q emission line has a complex multicom- 
ponent structure. Clearly, two or more systems of gaseous 
clouds, with different velocities along the line-of-sight, are 
present in this abnormal region. The non-circular motions 
apparently are not related to the bar dynamics because 
they are not symmetric relative to the galaxy center. We 
suppose an extragalactic gas cloud falling into the galactic 
disk. Similar shapes of the emission line profiles (at much 
larger distances from the center) were observed by us in 
NGC 1084 (Moiseev, 2000). 

4.3. NGC 2273 (Mrk 620) 

Van Driel & Buta (1993) have described a complex mor- 
phology of this galaxy: a large-scale bar of ~ 50" in diam- 
eter and two systems of spiral arms ending at distances 
l.'l and 1'.6 from the center. The influence of this large- 
scale bar appears clearly in the IFP velocity field in the 
H Q emission line (Fig. 14): the dynamical axis deviates 
by more than 20° from the line of nodes PA n at radii 
r = 3 — 45". The turn of PA dyn occurs in the oppo- 
site direction with respect to the PA of the isophotes at 
these distances (Fig. 2), in accordance with a common 
picture of gas flows in barred galaxies. The residual line- 
of-sight velocities along the minor axis have amplitudes 
of more than 50 — 80kms _1 (after pure circular model 
subtraction) at radial distances of ±10 — 13". The sign of 
the residual velocities corresponds to the gas inflow mo- 
tions 2 due to the large-scale bar. The PAdyn radial vari- 
ations in the H^ velocity field (from the MPFS observa- 
tions) are in agreement with those in the H Q velocity field. 
Nevertheless, the amplitude of the Hg rotation curve (at 
r = 1 — 6") is less by ~ 50kms _1 than the H Q rotation 
curve, and almost coincides with the rotation velocity of 
stars (Fig. 2). Apparently, the relatively weak H3 emission 
line is distorted by the broad absorption line from the un- 
derlying stellar spectrum. 

2 We assume that spiral arms are trailing. Thus the NW-side 
of the galaxy is the nearest to an observer. 



Moiseev et al.: Structure and kinematics of candidate double-barred galaxies. 



9 



Fig. 4. The same as Fig. 1 for NGC 2950 (To be seen in landscape) 
Fig. 5. The same as Fig. 1 for NGC 3368 (To be seen in landscape) 
Fig. 6. The same as Fig. 1 for NGC 3786 (To be seen in landscape) 



The [OIII] velocity field is significantly different from 
the measurements in the Balmer lines. At r < 3" the 
[OIII] line-of-sight velocities are constant, therefore the 
PAdyn and V ro t are uncertain in this region. At larger dis- 
tances from the center the [OIII] gradient of the line-of- 
sight velocities is smaller than that of the hydrogen lines. 
We suggest that the existence of non-circular motions in 
forbidden lines is caused by the jet propagating outward 
from the Seyfert 2 nucleus. Ferruit et al. (2000) found a 
jet-like feature in the HST [OIII]-images and ionized maps 
([OIII]/H Q ). It is seen, up to ~ 2" to the east from the 
nucleus and coincides with the extended radiostructure. 
The spatial geometry of the gas motions around the jet 
cannot be restored from the MPFS data, most probably 
these motions are outside the galactic plane. 

The stellar velocities do not reveal any non-circular 
motions in contrast with the gaseous velocities. The ellip- 
tical structure, elongated in E-W direction, appears in the 
stellar velocity dispersion map with a small drop of by 
10 — 15kms _1 in the rotation center (Fig. 2). 

Mulchaey et al. (1997) have discussed a secondary bar 
at r < 8" from their NIR-photometry data. The high- 
resolution HST images showed that this "bar" breaks up 
into several elliptical arcs with radii of r = 3 — 4" (see the 
residual brightness image in Fig. 2). Ferruit et al. (2000) 
described these features as a "pseudo-ring", and we sug- 
gest that it is a symmetrical spiral twisted into a ring. 
This conclusion agrees with the opinion of Erwin & Sparke 
(2002, 2003). However Petitpas & Wilson (2002) returned 
toward a secondary bar assumption in the central kpc of 
NGC 2273 based on their radio observation. We will dis- 
cuss the molecular gas properties of the sample galaxies 
in Sect. 5.7. It will be shown that in the case of NGC 2273 
the CO morphology agrees with an assumption of the in- 
ner stellar-gaseous disk embedded into the large-scale bar. 
The disk lies in the galactic plane and possesses a circum- 
nuclear spiral, as is mentioned above. 

4.4. NGC 2681 

The sharp variations of the isophotal parameters (Fig. 3) 
reflect the complex morphological structure of the galaxy. 
Friedli et al. (1996) connected the radial turns of the 
isophotes with two bars, first proposed by Wozniak et al. 
(1995). Erwin & Sparke (1999) offered a model of the cir- 
cumnuclear structure of NGC 2681, consisting of three 
independent bars, with position angles of the major axes 
PA 1 = 30° (20" < r < 60"), PA 2 = 75° (4" < r < 20") 
and PA 3 = 20° (r < 4"). The smallest bar- like structure 
is seen in the HST NICMOS images (Fig. 3). 



However, only isophotal analysis results are not enough 
to prove triple-barred structure because the photometry 
of the galaxy is a puzzle. We must note two items. 

Firstly, this lenticular galaxy has no global spiral arms 
but its optical color maps (Wozniak et al., 1995; Erwin & 
Sparke, 1999) reveal short tightly wound spiral arms and 
numerous dust lanes in the region of the "middle bar" (r = 
10 — 20"). The majority of the dust lanes are located to 
the south of the nucleus and are extended along the spiral 
arms. These flocculcnt mini-spirals also can be seen in 
our model-subtracted V-band image (Fig. 21) and unsharp 
masked R-band image by Erwin & Sparke (1999). 

Secondly, the disk's position angle is uncertain from 
available imaging data. It varies from PAq = 50° (Friedli 
et al., 1996) to 130- 140° (Erwin & Sparke,1999; 2003). 
The outer V-band isophotes (Fig. 3) have PA = 110 — 
120°. Also the PAd yn in the circumnuclear region (Fig. 3) 
differs from all these values. We tried to use a Fourier ex- 
pansion method to search for disk orientation parameters 
and failed completely, because there are no symmetrical 
spirals in the galaxy. 

The long-slit data help to determine the disk orienta- 
tion. Fig. 19 shows that the line-of-sight velocity gradient 
is absent along PA = 61°. Hence this is the direction of the 
minor axis (here we neglect possible non-circular motions 
in the bars). For more precise analysis we fitted velocities 
in all cross-sections by a model of pure circular rotation. 
The following values were found from a ^-minimization: 
PA = 148 ±5° and i = 25 ±5°. These parameters roughly 
agreed with the orientation of the outermost isophotes 
PA = 140° and i = 18° (Erwin & Sparke, 2003, they 
do not quote an uncertainty, but it is probably at least 
3°). 

The panoramic spectroscopy data (Fig. 3) show that 
stars at r > 3 — 4" exhibit a circular rotation where the 
value of PAdyn matches the disk's line of nodes. But at 
r < 2" we observe strong changes of the dynamical axis of 
the stars by 30°— 90°. The PAdyn variations in the velocity 
field of the ionized gas are smaller, but differ systemati- 
cally from that of the stars by 10—15° at r > 3". Also the 
symmetry center of the gas velocity field is shifted by 1.5" 
to the east from the continuum center, similarly to NGC 
470. In addition, a sharp turn of the isovelocities is seen 
at r < 2" (80 pc). So this region is dynamically decoupled 
concerning both the stars and the gas. At these distances, 
the value of PAd yn of the stars is closer to the orientation 
of the HST inner isophotes (PA = 25 - 35°). Such a re- 
lation between the "photometric" and "kinematic" axes 
does not point in the direction of the bar hypothesis, but 
of disk-like kinematics. Thus an inclined (or polar) stellar- 



10 



Moiseev et al.: Structure and kinematics of candidate double-barred galaxies. 



Fig. 7. The same as Fig. 1 for NGC 3945 (To be seen in landscape) 
Fig. 8. The same as Fig. 1 for NGC 4736 (To be seen in landscape) 
Fig. 9. The same as Fig. 1 for NGC 5566 (To be seen in landscape) 



gaseous disk exists in the central few tens parsecs of NGC 
2681. This disk is seen on HST- images of the galaxy from 
UV (WFPC2) to NIR (NICMOS). Erwin & Sparke (1999) 
accepted this disk as a third (inner) bar on HST-images. 
Unfortunately, the spatial sampling of the MPFS data is 
not sufficient for more certain conclusions on kinematics 
in this region (the seeing was 1.7"). 

Cappellari et al. (2001) have proved that NGC 2681 
hosts a nuclear starburst of <~ 10 9 yrs old. They consid- 
ered a merger process with a small gaseous companion as 
the most probable trigger of the nuclear star formation. If 
the companion had a spin direction inclined to the galac- 
tic plane, then the remnants of the merged gas will be ob- 
served as a central inclined disk, and the non-orthogonal 
disk will precess toward the galactic plane. However, a 
polar disk may be stable enough and rather long-living 
(Arnaboldi & Sparke, 1994). From this point of view the 
central structure of NGC 2681 is similar to that of galaxies 
with inner polar rings, such as IC 1689 (Hagen-Thorn & 
Reshetnikov, 1997) or UGC 5600 (Shalyapina et al., 2002). 
Nevertheless, the polar disk in NGC 2681 is smaller, of 
about 100 — 300 pc in diameter. 

The merger hypothesis makes it possible to explain 
the spiral arms and dust lanes at r < 20", which have 
a shape differing from ordinary dust lanes in large-scale 
spiral arms. Also these flocculent non-symmetrical spirals 
differ from the symmetrical two-armed dust spirals which 
are usually observed near a shock front in barred galax- 
ies. The H a +[NII] ionized gas emission lines is also ob- 
served only at r < 15" (Wozniak et al., 1995). Thus, the 
stellar-gaseous disk, with spiral arms and perturbed dust 
layer 3 , is nested in the center of the large scale-bar at 
r < 15 - 20" (600 - 900 pc). A strong asymmetry in the 
dust distribution may be associated with a possible warp 
of the gaseous-dust disk relative to the stellar one. We 
suggest that all the features mentioned above (inner po- 
lar disk, abnormal shape of the spirals, and warp of dust 
layer) are caused by a recent merging event. 

We try to estimate the positions of ILRs in NGC 
2681. Figure 20 shows the polynomial approximation of 
the stellar rotation curve. However, the amplitude of cir- 
cular velocities must be larger than these streaming veloc- 
ities because velocity dispersion provides additional sup- 
port against gravity. For asymmetric drift correction we 
used Eq. (2) from Aguerri et al. (2003) and their assump- 

3 We suggest that the spiral in NGC 2681 probably has a 
dust origin. Indeed, the spiral appears in the optical images 
and is absent in the J-band residual image (Fig. 3) as well as 
in NIR color maps (Friedli et al., 1996) 



tions (the ratio of azimuthal to radial components of the 
velocity dispersion o^/ctr = 1/^2), <t* and a decrease of 
the volume density of the disk that is exponential with ra- 
dius). The drift correction was calculated for az I or values 
0.7 and 1.0, bracketing the possible range for early-type 
galaxies (Aguerri et al., 2003). The corrected circular ve- 
locities and the resonance curves shown in Fig. 20 are 
truncated at the radius of 15" because the asymmetric 
drift approximation fails in the bulge-dominated region 
(see surface brightness profile in Fig. 21). For the calcu- 
lation the angular velocity fip of the bar we make the 
following assumptions: the maximal circular velocity is 
about 220 - 260 km s" 1 ; the bar size is R b = 60"; the 
ratio of the radii can be 1 to 1.6 Rb, this is the range ob- 
served for bars in lenticular galaxies (Aguerri et al., 2003). 
The probable range of the large-scale bar angular velocity, 
fip = 35 — 70 kins" 1 kpc -1 , is marked as gray rectangular 
in Fig. 20. This figure also demonstrates that one or two 
ILRs could exist near the radius r w 20", if the fip is 
smaller than the value mention above. This distance coin- 
cides with the approximate size of the inner dust/gaseous 
disk described above. In this case the gas was accumulated 
by the gravitational forces of the bar. Of course we must 
remember that the resonance curves are calculated for a 
symmetric potential without any bar perturbations. But 
we think that our evaluation of the ILR radius can be used 
as a zero approximation of the position of the resonances. 

In conclusion, we suggest the following interpretation 
of the global morphology of NGC 2681. The large-scale 
low-contrast bar dominates in the optical and NIR im- 
ages at r = 50 — 60" and has the ILRs at r sa 15 — 25". 
The decoupled gaseous disk is observed inside the reso- 
nance region. This inner disk is elliptical (non-circular) 
because the elongated orbits of stars change their major 
axis orientation at the resonances (Lindblad, 1999). This 
elliptical disk was accepted by Erwin & Sparke (1999) as 
a secondary ("middle") bar. The recent dwarf compan- 
ion merging leads to the structure perturbation of the 
kiloparsec-size disk with inclined (polar) orbits in the cen- 
tral 100 pc. 

4.5. NGC 2950 

Our value of PAq differs slightly (by 10°) from the mea- 
surements of Friedli et al. (1996). The radial dependences 
of the isophotal PA and e, in all observed bands, coin- 
cide with those of Wozniak et al. (1995) and Friedli et al. 
(1996). The large-scale bar (r « 30 - 40" in size) is seen 
over the distributions of the m = 2 Fourier harmonic and 
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Fig. 10. The same as Fig. 1 for NGC 5850 (To be seen in landscape) 



Fig. 11. The same as Fig. 1 for NGC 5905 (To be seen in landscape) 



Fig. 12. The same as Fig. 1 for NGC 6951 (To be seen in landscape) 
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Fig. 19. NGC 2681: Stellar velocities measured along different 
position angles. The solid line marks the best-fit model of the 
rotation curve 



has a turn of isophotal PA accompanied by an ellipticity 
peak (Fig. 4). 

Wozniak et al. (1995) and Friedli et al. (1996) ex- 
plained the inner PA turn, and the peak, at r < 6", as 
a secondary bar. Surprisingly, this "bar" does not distort 
the circular rotation of stars in the circumnuclear region. 
The PAdy n of the stellar velocity field coincides with PA 
of the line of nodes. 

We do not detect any non-circular motions in this ve- 
locity field despite the strong turn of isophotes at these 
radii. There are sharp elliptical structures decoupled in 
the stellar velocity dispersion map at r = 1 — 5" (Fig. 4). 
It is extended along PA w 150°, which is the direction of 
the major axis of the primary (external) bar. 



4.6. NGC 3368 (Messier 96) 

Sil'chenko et al. (2003) analyzed the results of panoramic 
spectroscopy of this galaxy including our observational 
data. However, that work is mostly focused on the stellar 
population properties. We will briefly describe the main 
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Fig. 20. top Rotation velocities of the stars in NGC 2681. 
Different symbols show MPFS (filled) and long-slit (open) 
data. The solid line is a polynomial fitting. The other lines 
mark the velocity streaming curves where asymmetric drift 
terms are taken in account with <t z /(Tr = 0.7 (dotted) and 1.0 
(dashed), bottom The angular velocity (fi) and il—n/2 (where 
k is the epicyclic frequency) curves calculated from the differ- 
ent models of the rotation curve. The thick rectangle shows 
the possible range of the angular velocity of the bar. 



properties of the inner kinematics of this galaxy, which is 
a member of the Leo I group. 

Observations with the scanning IFP allow us to map 
the large-scale velocity fields of the gas in two emission 
lines (Fig. 15). The H Q emission is mostly concentrated in 
the starforming ring, at distances of 50 — 70" from the cen- 
ter. Unfortunately, the stellar absorption features and the 
overlapping interference orders distort the emission line 
profiles in the H Q data-cube, over the whole central region 
of the galaxy. We cannot resolve this problem because the 
free spectral range of the IFP is rather small (28 A). In 
contrast, the [Nil] velocity field has been constructed over 
the full range of radii. The dynamical position angle of 
the outer gaseous disk, PA dyn = 170°, at r = 30 — 200", 
differs strongly from the orientation of the line of nodes of 
the stellar disk, PA = 135°. The last value was calculated 
from the Fourier expansion of the NIR images ( "spiral cri- 
terion" , Sect. 3). So the line of the maximum of the m = 2 
harmonic has a regular spiral shape and fits both spiral 
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Fig. 13. The same as Fig. 1 for NGC 7743 (To be seen in landscape) 

Fig. 14. IFP observations in the H a line, left - NGC 470. At the top continuum and H a images plotted on a square-root scale. 
The cross marks the dynamical center. At the bottom the velocity field with isovelocities is shown. The thick black contour 
corresponds to the systemic velocity, while the thin contours indicate line-of-sight velocities in steps of ±50 kms -1 . The last 
plot shows the rotation curve and the kinematic position angles. Here the dotted line marks the PAo value. There is a ghost 
from a bright star located in position (20", 35"). right - same for NGC 2273. (To be seen in landscape) 

Fig. 15. NGC 3368. IFP observations in the [Nil] (left) and H Q lines (right). The dotted line correspond to the Position Angle 
of the gaseous disk of 170° (see text). (To be seen in landscape) 



arms of the residual brightness distribution (Fig. 5) only 
for this value of PA . The isophotes, distorted by these 
spiral arms, were erroneously interpreted by Jungwiert et 
al. (1997) as a large-scale primary bar. 

The orientation of the ionized gas disk PAd yn = 
170 — 175° is in good agreement with the outer part of 
the neutral hydrogen velocity field (Schneider, 1989) and 
with CO kinematics at r < 15" (see Fig. 1 in Sakamoto et 
al., 1999). Thus, we suggest that the gas (ionized/neutral) 
and the stars in this galaxy rotate in different planes. All 
the gas in this galaxy, which is the brightest member of the 
western half of Leo group, has an external origin, in agree- 
ment with Schneider (1989). See Sil'chenko et al. (2003) 
for detailed arguments in favor of this possibility. 

Good agreement exists between the measurements 
of the PAdyn from two independent IFP velocity fields. 
Nevertheless, the [Nil] rotation velocities are systemat- 
ically lower, by ~ SOkms" 1 , than the H^ones. If the 
gaseous disk is inclined to the main galaxy plane, then 
shock-wave fronts can develope at the cross-section of the 
global stellar and gaseous disks, because the gas strikes the 
gravitational well. Thus, the low [Nil] velocities may be 
explained if the shock-excited gas, which is slowed down 
by collisions with the stellar disk and spiral arms (these 
intensify the contrast of the gravitation potential), emits 
mostly in the forbidden emission line. The discrepancy be- 
tween the velocities in the different emission lines is not an 
artificial instrumental effect. NGC 4736 was also observed, 
with the same technique, in H Q and [Nil] lines during the 
same observation run. In the NGC 4736 data cubes the 
velocities in both emission lines are consistent with each 
other. 

There is no large-scale ( "primary" ) bar in this galaxy, 
as it was mentioned above. However, a small-scale (r w 
5") nuclear mini-bar is present and extends along PA = 
120° — 125°. This bar provokes a turn of PAd yn in the 
stellar velocity field by ~ 25°, from the line of nodes in 
the opposite direction with respect to the PA turn of the 
inner isophotes (Fig. 5). The central asymmetry of the cr» 
distribution is confirmed by long-slit spectroscopy of Vega 
Beltrafi et al. (2001). 

The sharp turn of the isophotes, in the HST images 
inside r = 2", may be associated with strong dust lanes. 
Each of these lanes crosses the images near the nucleus and 
can be seen both in the optical and NIR images (Fig. 5). In 



Sil'chenko et al. (2003) we argue that a gaseous-dust mini- 
disk is located in the polar plane of the bar. Unfortunately, 
our MPFS and IFP data have too low angular resolution 
to study the polar disk kinematics. 

4.7. NGC 3786 (Mrk 744) 

Our disk orientation parameters for this galaxy are in 
agreement with the results of Afanasiev & Shapovalova 
(1981). Afanasiev et al. (1998) from their imaging data 
found two bars with semi-major axes of 7" and 25" in 
this Sy 1.8 galaxy. However, our Fourier analysis of the 
optical and NIR images shows that the line of the max- 
imum of the m = 2 harmonic, at r > 6 — 7", is wound 
into a regular spiral that coincides with the global two- 
armed spiral pattern of the galactic disk. This spiral also 
appears in the residual brightness maps (Fig. 6). Thus, 
the radial variations of PA of the inner isophotes caused 
by the spiral arms have led Afanasiev et al. (1998) to an 
erroneous conclusion about a large-scale bar in NGC 3786. 
The situation is the same as in the case of NGC 3368. 

The PAdyn of the velocity field of stars at r < 6", 
differs from the PAq of the line of nodes by more than 10°. 
The deviations occur in the opposite direction relative to 
the central isophotes (Fig. 6). This indicates a dynamically 
decoupled central mini-bar of 2 kpc in diameter. The line 
of the maximum of the m = 2 harmonic at r < 5 — 6" has 
a constant azimuthal angle, a further confirmation of the 
existence of the mini-bar. The mini-bar was predicted by 
Afanasiev & Shapovalova (1981) from their long-slit data 
and confirmed photometrically by Afanasiev et al. (1998). 

The ionized-gas velocity fields in the Hg and [OIII] 
lines arc different. More precisely, the orientations of 
PA 

dyn in velocity field of the stars and gas in H^ coincide, 
but the position angle, calculated from the [OIII] velocity 
field of the gas, deviates from the stellar one. The devia- 
tions may be due to the presence of shock fronts at the bar 
edges, formed in the disk gas under the action of the grav- 
itational well of the bar. The post-shock gas decelerates 
and emits in forbidden lines (Afanasiev & Shapovalova 
1981, 1996). The rotation velocities in the [OIII] line are 
smaller than those in the Hp (Fig. 6). This feature may 
also be connected to braking shock-excited gas at the bar 
edges. 
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Fig. 16. IFP observations of NGC 3945 in the [Nil] line. 

Fig. 17. NGC 4736. IFP observations in the [Nil] (left) and H a lines (right). (To be seen in landscape) 
Fig. 18. NGC 6951. IFP observations in the [Nil] (left) and H a lines (right). (To be seen in landscape) 



The stellar velocity dispersion field has drops by ~ 
50kms _1 in the center compared with ct* at r — 3 — 4". 
This feature will be discussed in Sect. 5.6. 

4.8. NGC 3945 

The radial variations of the isophotal parameters over all 
observed bands are in accordance with early observations 
of Friedli et al. (1996) and Erwin & Sparke (1999). The 
global spiral pattern is absent in this galaxy. The param- 
eters of the disk orientation were calculated under the as- 
sumption of a circular shape of the outer low-brightness 
ring (r = 120 — 150"). The values of PAq and io agree 
(within the errors given in Table 7) with the results of 
Erwin & Sparke (1999) and Erwin et al. (2003). 

A large-scale bar extends along the galactic minor axis. 
Erwin & Sparke (1999) accepted that the circumnuclear 
disk nests in the bar, because the isophotal PA matches 
PAq at r = 4 — 18". They also favored the existence of 
a secondary small bar embedded in the disk, and found a 
low-contrast ring at r = 7" as well. This ring is seen in 
our residual HST- images (Fig. 7). The origin of this ring 
is discussed by Erwin et al. (2001). 

The velocity field of the stars in the circumnuclear disk 
shows a regular circular rotation (PAd yn coincides with 
the PAo). That the kinematics is disk- like agrees with the 
conclusion by Erwin et al. (2003) that the inner disk is 
more luminous (and massive) than the bulge in this region. 
Also we find no any non-circular motions of the stars due 
to a possible secondary bar at r < 3" suggested by Erwin 
& Sparke (1999). The rotation curve agrees with early 
measurements of Bcrtola et al. (1995). A small difference 
in the amplitude of the rotation curves, observed in 2000 
and 2002, is probably due to a velocity gradient, smoothed 
by beaming effects, in observations with a worse seeing in 
2000 (Table 1). 

The velocity field of the ionized gas is puzzle. Bertola 
et al. (1995) found [OII]A3727A line emission only in the 
very central region (r < 3") and showed that the gas co- 
rotates with the stars. Nevertheless, our MPFS data con- 
tradict their conclusions. First of all, we confidently detect 
the [Nil] line emission up to a distance of 10" from the cen- 
ter. Secondly, inside the radius of 6" (about 0.5 kpc), the 
line-of-sight velocities of the gas are inverted with respect 
to the stellar velocities. Thus, a counter-rotating gaseous 
disk is located at r < 5 — 6". 

The direction of gas rotation becomes coterminous 
with the stellar rotation at larger distances from the 
nucleus (sec Fig. 7). The large-scale IFP velocity field 
(Fig. 16) confirms the fact of normal gas rotation at large 



radii, up to 140" (~ 11 kpc). The variations of PAd yn a,t 
r = 5" — 35" in the IFP velocity field are not real. This 
artifacts are caused by spectral overlapping of the strong 
stellar H Q absorption and low-contrast [Nil] emission lines 
in the central region. The small spectral range of the IFP 
observations does not allow us to separate absorption and 
emission lines. In contrast, at r = 115" — 140" the stel- 
lar continuum has a low contrast, while the [Nil] lines are 
bright since they are emitted by some HH-regions in two 
spirals twisted into the outer ring. 

4.9. NGC 4736 (Messier 94) 

The outermost isophotes of our NIR images are limited to 
r w 50", and we cannot calculate the parameters of the 
outer disk orientation from these images. The application 
of the "tilted-ring" model to the [Nil] velocity field at 
r = 80 — 120" makes it possible to obtain the PA value. 
The disk inclination is small and uncertain from velocity 
field analysis. Mulder (1995) found io = 31° from the outer 
optical isophotes. We accept this value throughout our 
paper. 

The stellar velocity field shows a regular circular ro- 
tation with PA dyn = 300 - 305° (Fig. 8). These values 
differ slightly from the accepted PA , but fall inside the 
PAq errors (Table 7). From long-slit observations, Mulder 
(1995) obtained PA dyn = 285° - 290° for the stellar com- 
ponent in the galactic disk as a whole. This value differs 
certainly from our measurements at r < 10". 

Our H Q velocity field is similar to the data of Mulder 
(1995), but it has no measurements inside r < 25". Here 
the broad H Q absorption line distorts the emission line pro- 
file. The absorption does not affect the [Nil] emission. 
The [Nil] velocity field reveals non-circular ionized gas 
motions, as is shown from the fact that PAd yn turns at 
r = 5 — 30". The non-circular motions are most likely con- 
nected with bar- like structure at r < 15 — 20", detected 
from isophotal parameters variations: a turn of PA and 
a e peak (Fig. 8). This conclusion confirms the result by 
Wong & Blitz (2000) who suggest that molecular gas kine- 
matics is consistent with radial flows in the central stellar 
bar (~ 30" in diameter). 

Shaw et al. (1993) also considered a twist of NIR 
isophotes at r < 5 — 7", which may be explained as a sec- 
ondary bar effect. We must notice that the most highest- 
contrast part of the circumnuclear mini-spiral is located 
at these radii and distorts the shape of the isophotes in 
the HST images. This spiral structure is studied in recent 
work by Elmegreen et al. (2002). 
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The dust lanes in the spiral appear in the HST image 
up to r = 20 — 30". Thus, the deviations of PAd yn of the 
ionized gas velocity field from PA at these radii can be 
explained by non-circular gas motions in the mini-spirals 
instead of motions inside a large-scale bar as was men- 
tioned above. 



4.10. NGC5566 

The disk inclination is considerably smaller than the value 
of io = 80° measured by Jungwiert et al. (1997). We found 
io = 64° by means of the "spiral criterion" . This value 
is in good accordance with the outer isophote ellipticity 
(Fig. 9) and with Kent's (1984) CCD-observations. The 
large-scale bar of 25" in size is elongated at PA w 155°. 
Two strong spiral arms start from the bar ends. The ve- 
locity field of stars lacks significant non-circular motions 
as PAdyn coincides with PA within the errors. The stel- 
lar rotation velocity at r = 10 — 14" reaches 250 km s -1 , 
which corresponds to the HI rotation curve of Helou et al. 
(1987). 

Jungwiert et al. (1997) explained the isophotal ellip- 
ticity peak and the PA turn at r < 5" as secondary bar 
features. We disagree with this interpretation for the fol- 
lowing reasons. First of all, the position angles of the 
"secondary bar" and the line of nodes are very similar. 
Secondly, two symmetrically twisted spiral arms appear 
in the JHK 7 residual images. This distribution of the un- 
derlying stellar component is not characteristic for a bar 
(Fig. 9). And finally, a grand-design two-armed nuclear 
spiral is seen also in the HST residual image (Fig. 9) at 
these radii. 

We must note that a slowly rotating mini-bar can pro- 
duce a spiral-like response in the disk outside the ends of 
the bar, where an ILR must be located, as has been shown 
in the numerical simulations by Combes & Gcrin (1995) 
and in the analytic estimations by Fridman & Khoruzhii 
(2000). But in the case of NGC 5566 the mini-spiral devel- 
ops inside the radius, where a possible nuclear mini-bar 
is proposed by Jungwiert et al. (1997). Secondly, we find 
no non-circular motions produced by a possible nuclear 
bar in the stellar velocity field, because PAd yn agrees with 
the disk PAq (see above). And thirdly, a possible bar-like 
structure is not visible in the residual brightnesses depro- 
jected to the galactic plane (see Fig. 21). 

For all these reasons we suggest that the mini-spiral is 
not connected with a possible secondary bar. Therefore, 
NGC 5566 is an ordinary barred galaxy. It has a circum- 
nuclear disk (r = 5" — 8", or 500 — 750 pc in linear radius) 
with a mini-spiral nesting inside this bar. The properties of 
the inner isophotes (PA — 205 — 210°, e max = 0.55) indi- 
cate that this kiloparsec-size inner disk lies in the galactic 
plane. The mini-spiral is traced by dust lanes but the ma- 
jority part of the dust clouds is concentrated outside the 
spiral arms in the NW half of HST PC2 image. 



4.11. NGC 5850 

Since our NIR images have a field of view that is not 
large enough to cover the entire galaxy, we used the R- 
band image, from the digitized atlas of Frei et al. (1996) 
to determine the orientation parameters. The spiral arms 
are low-contrast outside the ring around a large-scale bar 
(r > 80"). Higdon et al. (1998) have shown that the outer 
spirals have a mixed (2- and 3-armed) structure. However, 
the orientation parameters, determined by means of the 
"spiral criterion", from the m — 2 harmonic are in good 
agreement with preliminary data of Friedli et al. (1996), 
and with the kinematics of neutral hydrogen if the possible 
warp of the Hi-disk is taking into account (Higdon et al. , 
1998). Since the inner isophotes have a sharp turn at r < 
10", Friedli et al. (1996) have assumed the existence of a 
secondary bar in this galaxy. 

The MPFS data show that PA dyn of the stellar rota- 
tion coincides with PA , while in the ionized gas PAd yn , 
in the circumnuclear region, differs by more than 50 — 60° 
from the line of nodes position angle (Fig. 10). Note that 
the ionized gas PAd yn almost aligns with the PA of the 
inner isophotes. This radial behavior of PAd yn is typical 
for a disk inclined to the galactic plane. Nevertheless, if 
we assume that the gas motions are produced in the galac- 
tic plane, they would correspond to a radial outflow from 
the nucleus 4 with velocities of 50 — 70 kms -1 , without as- 
suming any circular rotation. Such outflows are typical 
of Seyfert galaxies, nevertheless, no AGN or starburst- 
like features are observed in the optical spectrum of the 
galaxy nucleus (see discussion in Higdon et al., 1998). A 
more reasonable assumption is that the gas, at r < 6 — 7", 
rotates in a polar plane with respect to the global galactic 
disk. In this case, the polar gaseous disk lies roughly in 
the smallest principal-plane cross-section of the outer bar, 
orthogonal to the major axis of the bar, as in NGC 3368 
(Sect. 4.6). The hypothesis of a polar disk is supported 
by the fact that NGC 5850 has undergone a recent col- 
lision with the nearby galaxy NGC 5846 (Higdon et al. 
1998). Through their interaction, part of the gas could be 
transported to polar orbits. What is the inclination of this 
disk to the line-of-sight? It is easy to show that the angle 
between the inner disk plane and the outer galactic disk, 
Ai, satisfies the following relation: 

cos Ai — ±sinii sini cos(PAi — PA n ) + cos?i cosio, (2) 

where i\ is the angle between the line-of-sight and the 
plane perpendicular to the disk, and PA\ is the posi- 
tion angle of the major axis of the disk. By substituting 
PA a = 335°, i = 37° (Table 7), and PA t = PA dyn = 35° 
into this expression, we predict that the disk will be polar 
(Ai = 90 ± 5°) for h = 65° - 74°. The gas rotation curve 
in Fig. 10 was calculated for %\ = 65°. 

At r > 5", the gas rotation velocities may be lower 
than the stellar ones, but both velocities are consistent 

4 Here, we use Higdon's et al. (1998) suggestion about the 
galaxy disk orientation, namely that its western half is closest 
to the observer. 
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within the errors. An extrapolation of the HI rotation 
curve (Higdon et al., 1998) to the inner radii gives V rot ~ 
30 — 45 kms -1 at r = 5" — 8", in good agreement with our 
ionized gas measurements. 

4.12. NGC5905 

Unfortunately, we have not obtained the image of the 
outer isophotes, extended beyond a large-scale bar at 
r w 37" (the bar size is taken from Friedli et al., 1996). 
Therefore, we accepted io — 40°, which was derived by 
van Moorsel (1982) from the HI kinematics analysis. Also 
we calculated PA n = (130 ± 5)° from inspection of the 
HI isovelocities map (Fig. 8 in van Moorsel, 1982). The 
last value agrees perfectly with the results of deep CCD- 
photometry by Kent (1984). 

The Hp velocity fields of stars and ionized gas exhibit 
PAdyn deviating from PAq at r = 1" — 9", in the opposite 
direction with respect to the PAs of the inner isophotes. 
Hence they should be caused by non-circular motions in 
the large-scale bar. The observed slower [OIII] ionized gas 
rotation compared with the Hp value and the difference of 
the PAdyn values in these lines are most likely connected 
with shock fronts in the bar edges, as in the case of NGC 
3786 (Sect. 4.7). 

Friedli et al. (1996) showed that the turn of the in- 
ner isophotes, at r < 10", is caused by a secondary 
bar. Nevertheless, they used an obviously wrong value of 
PAq = 45°. Note that a relatively small value of the el- 
lipticity of the inner isophotes, and the coincidence of the 
isophotal PA with PAq does not allow us to explain the 
central structure with a secondary bar. Most probably, the 
twist of the central isophotes at r = 5" — 10" is caused by 
the transition of the surface brightness distribution char- 
acter from a bar-dominated to a central prolate bulge or 
circumnuclear disk. This situation is similar to that of 
NGC 5566 (Sect. 4.10), but without introducing the mini- 
spirals. 

4.13. NGC 6951 

The orientation parameters from Table 7 agree with 
Perez's et al. (2000) results. The ionized gas velocities, 
measured by them in three different positions of the slit, 
coincide completely with our velocity fields. However, we 
did not detect the two kinematic components in stellar ro- 
tation, perhaps because of the low MPFS spectral resolu- 
tion, relative to that of Perez et al. (2000) long-slit data. In 
the star velocity field, the PAd yn differs, by 5° — 10°, from 
the PAo (Fig. 12). The turn of the dynamical axis is not 
very significant. This fact is connected with non-circular 
motions in the external bar, because the bar isophotes 
deviate in the opposite direction relative to PA . 

At large distances from the center, the IFP H a and 
[Nil] velocity fields are similar (Fig. 18). Nevertheless, 
the [Nil] rotation velocities are lower than those of the 
H Q inside the circumnuclear starforming ring at r = 



5" - 10", while the IFP H Q and [Nil] and the MPFS 
H^ PAdyn variations are the same. In contrast, the MPFS 
[OIII] velocity field exhibits more significant non-circular 
ionized gas motions. The dynamical axis, found from this 
line measurements, turns by ~ 30° at r < 3", and the ro- 
tation velocity is lower than the Hp velocities by a factor of 
2. The intensity of the [OIII] emission decreases sharply 
outside this region. It is possible that radial motions of 
highly excited ionized gas, with a speed of the same order 
of the rotation velocities, occur here. Radial motions are 
connected with the outflow from an active nucleus, which 
is of a transition type between a LINER and a Sy 2 (Perez 
et al., 2000). Moreover, the non-thermal radio continuum 
image shows a jet-like feature, extended from the nucleus 
up to 80 pc (Saikia et al., 2002). Unfortunately, the MPFS 
data cannot resolve this asymmetric structure. 

There is a circumnuclear ring at r = 5 — 7", nested 
in the large-scale bar, which has a r « 60". The turn 
of isophotes, in the ground-based optical images, inside 
r < 6" was explained by Wozniak et al. (1995) with a pos- 
sible secondary bar. However, Friedli et al. (1996), based 
on NIR-photometry, have assumed that the isophote dis- 
tortion is caused by a complex distribution of dust and 
star-forming regions inside the central kiloparsec. The 
HST image resolved this ring into individual starforming 
knots. The ring seems to be elliptical, but it is almost cir- 
cular after deprojection onto the galactic plane (PA and 
e of the central isophotes correspond to the outer disk 
orientation). Indeed, the HST high- resolution images con- 
firm the lack of a secondary bar as an elongated stellar 
structure (Fig. 12). 

The HST residual brightness maps reveal a multi- 
armed flocculent spiral (Fig. 12), described in detail by 
Perez ct al. (2000) . The mini-spiral is probably associated 
only with the distribution of gas and dust rather than 
with the stellar component. As has been shown by Perez 
et al. (2000), the nuclear spiral structure, clearly seen in 
the V band, completely disappears in the H band, where 
the effect of dust absorption is much weaker. 

Figs. 12 and 18 reveal that the position angle of the 
dynamical axis, calculated from the H a , H^ , and [Nil] 
velocity fields, at r = — 8" deviates by 10 — 15° from 
the PAq. This proves the significant role of non-circular 
motions in the disk kinematics. Since there is no inner bar 
in this region, we can offer the following interpretation of 
the observed picture. A gaseous dusty disk of 6"— 8" (400— 
600 pc) in radius is embedded into the large-scale bar. In 
this disk, a multiarmed spiral structure has developed and 
perturbs the circumnuclear gas rotation. The dynamical 
decoupling of this disk is also confirmed by a high central 
molecular-gas density (Kohno et al., 1999), and by the 
location of two ILRs of the large-scale bar (Perez et al., 
2000). 

Based on their IFP observations, Rozas et al. (2002) 
presented similar arguments to explain the kinematic 
decoupling of the inner ionized-gas disk in NGC 6951. 
However, they tried to interpret the gas non-circular mo- 
tions in terms of an inclined disk or a dissipation of the 
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secondary bar. In a forthcoming paper we are going to 
consider a detailed kinematics and circumnuclear struc- 
ture of this galaxy, based on new 2D-spectroscopy data, 
with sub-arcsecond angular resolution (Moiseev, private 
communication) . 

4.14. NGC 7743 

A regular two-armed spiral is seen in the ground-based 
images, most clearly in residual images (Fig. 13). The turn 
of the isophotes, inside r w 50", was explained by Wozniak 
et al. (1995) as the combined influence of a bar and a 
triaxial bulge, with the bulge dominating at r < 10". Our 
analysis of residual images and the behavior of the m = 2 
Fourier harmonic show that the bar has a smaller size, and 
that the spiral structure distorts the elliptical isophotes 
at r > 20" - 30". The sharp isophote twist, at r < 3", 
can be seen as well in ground-based as in HST images 
(see Fig. 13). The WFPC2 residual optical images reveal 
a complex multi-armed spiral (Fig. 13). 

This spiral completely disappears in the near-infrared 
NICMOS images. Apparently, this spiral has a dust 
origin (like the gaseous-dusty spiral in NGC 6951). 
Unfortunately, our MPFS ionized-gas velocity field, in the 
mini-spiral region, is uncertain. The [OIII] emission line is 
weak and noisy. We can only conclude that at r < 2" — 4" 
the gas follows the stellar rotation. A region with non- 
circular motions is located at large distances, South from 
the center . 

The velocity field of the stars shows a regular rotation 
with almost constant PAd yn , while the mean PAd yn value 
differs, by <~ 30°, from PAq. The deviation goes in 
the same direction as the turn of the inner isophotes. 
Moreover, the PAdyn — 300 — 310° coincides with the ori- 
entation of the inner NICMOS isophotes. This is normal 
for an inclined rotating nuclear disk. Nevertheless, on the 
our ground-based and HST images this region (r < 8") 
is not distinguish, except by the mini-spiral at r < 4". 
It is possible that the accepted value of PAq is not cor- 
rect. Additional observations of this galaxy focused on the 
large-scale gas and stellar kinematics are needed. 

The stellar velocity dispersion field is noisy (Fig. 13), 
but it is in good agreement with the long-slit observations 
of Kormendy (1982), concerning the values of the central 
velocity dispersion, and the position of the tr* peak, shifted 
by 3" — 4" from the photometric center. 

5. Discussion 

5.1. Circumnuclear morphology of the candidate 
double-barred galaxies 

In Table 8 we collect all structural and kinematic features 
which have been found in the observed galaxies. The table 
contains galactic name, type (according to the RC3 cat- 
alogue), accepted distance from the HyperLeda database 
(H = 75kms _1 Mpc -1 ), bar major-axis length and posi- 



tion angle 5 and brief description of the structures at differ- 
ent radial scales. The properties for a single bar are given 
in each case because we could not confirm kinematically 
secondary bars in any galaxy. Usually the bars described 
were denoted as "primary" in the references. However in 
NGC 3368 and NGC 3786 only nuclear dynamically decou- 
pled bars are present (denoted as "secondary" in the refer- 
ences), and "large-scale primary bars" are, in fact, a con- 
sequence of isophotes distorted by spiral arms (Sect. 4.6 
and 4.6). 

We have detected various types of deviations of the 
stellar and gaseous motions from the normal regular rota- 
tion in the circumnuclear regions of the galaxies observed. 
These types are separated below. 

5.2. Non-circular motions in bars 

The PAdyn calculated over the ionized gas velocity fields 
disagrees with the PAq of the lines of nodes in all galaxies 
where emission lines have been detected. In Sect. 4 we have 
shown that the central turns of ionized-gas PAd yn occur in 
opposition to the "primary bar" isophotes independent of 
the orientation of the secondary bar-like structure (exclud- 
ing galaxies with inner polar and counter-rotated gaseous 
disks, Sect. 5.4 and 5.5). We explain such behaviour of 
the "dynamical axes" as non-circular gas motions (radial 
flows) in a large-scale bar (or nuclear mini-bars in the case 
of NGC 3368 and NGC 3786) in agreement with models 
described in Sect. 3. In the stellar velocity fields the am- 
plitudes of PAdyn deviations are usually less than those 
of the gas; in half of the sample PAd yn is close to the 
PA„ value. In NGC 3368, NGC 3786, NGC 5905, and 
NGC 6951 significant (by 7 - 20°) PA dyn deviations are 
detected. Again comparing these with the isophotal PA 
they are typical signatures of non-circular motions of the 
stars in a single-barred galaxy, or of radial motions due to 
elongated stellar orbits, including the case of twisted x 2 
orbits. 

NGC 2950 is an interesting case, where the inner 
isophotes are turned by more than 50° (Fig. 4), but the 
PAdyn for the stellar component is aligned with the line of 
nodes. Neither the outer bar nor the inner bar-like struc- 
ture have an appreciable effect on the stellar velocity field. 
It is possible that both disk stellar motions within the 
bar and the rotation of the bulge, whose contribution in 
this SB0 galaxy must be significant, are observed along 
the line-of-sight. Since the spectral resolution was too low 
to separate two dynamical components within the line- 
of-sight velocity distribution (LOSVD), the mean velocity 
field corresponds to circular rotation. The stellar compo- 
nent associated with the outer bar in NGC 2950 affects 
the velocity dispersion distribution, whose increase (ellip- 



5 The values of a and PA were estimated from our isophotal 
analysis using same method as in Wozniak et al.(1995). Erwin 
(2004) has provided new definitions for bar semi-major axis 
measurement. 
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Name Type D 1 Bar Circumnuclear structure 





(RC3) 


Mpc 


a," 


PA,° 


r 


It 


r, kpc 


type of the feature 


NGC 470 


SAb 


32.6 


30 


20 





- 6 


0- 


- 1.1 


lopsided nuclear disk 












8- 


-20 


1.3 


-3.2 


gas radial motions in the bar 


NGC 2273 


SBa 


25.7 


25 


108 





-3 


0- 


-0.4 


jet outflow in [OIII] line 












1 


- 4 


0.1 


-0.5 


nuclear disk with mini-spirals 












3 . 


- 27 


0.4 


-3.4 


o"fl^ rflrlifil TnrrHnnQ in flip liar 


NGC 2681 


SABO/a 


9.6 


60 


30 





- 5 


0- 


- 0.2 


stellar/gaseous polar (inclined) disk 












5- 


- 20 


0.2 


-0.9 


decoupled disk with dusty mini-spiral 


NGC 2950 


SBO 


18.9 


40 


155 












NGC 3368 


SABab 


10.6 


5 


125 





- 2 


0.1 


gaseous polar disk 















- 5 


- 


- 0.3 


stellar radial motions in the mini-bar 












5 - 


- 200 


0.3 


- 11 


global gaseous disk inclined to the stellar one 


NGC 3786 


SABa 


36.1 


6 


70 


2 


-6 


0.4 


- 1.0 


gas/stellar radial motions in the mini-bar 


NGC 3945 


SBO 


17.7 


40 


72 





-6 


0- 


- 0.5 


counter-rotating gaseous disk 


NGC 4736 


SAab 


4.8 


18 


20 


0- 


- 25 


0- 


- 0.6 


mini-spiral structure with non-circular gas motions 


NOC 5566 


SBab 


19.9 


20 


155 


o 


_ 8 


0- 


- 0.8 


miplpfiT" rli^k with Tnini-^nirfil 
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NGC 5850 


SBb 


34.2 


80 


115 





-6 


0- 


- 1.0 


polar gaseous disk 


NGC 5905 


SBb 


47.2 


30 


200 


o 


_ 8 


0- 


- 1.9 


Tiiiplpflr H i * pffiQ /Q+pllfir TflHial motions in tnp Vifir 

11 LH_ ±<_ Cll ttlLl/ O Lr^lltXl 1 C11_11C1j1 IIIUIIUIIlI ill LillC Ufll 


NGC 6951 


SABbc 


21.9 


80 


90 





-3 


- 


- 0.3 


jet outflow in [OIII] line 












1 


-8 


0.1 


-0.8 


stellar radial motions in the bar 












3 


-8 


0.3 


-0.8 


nuclear disk with mini-spiral and non-circular gas motions 












30 


- 55 


3.2 


- 5.9 


gas radial motions in the bar 


NGC 7743 


SBO 


24.4 


30? 


95 





-4 


0- 


- 0.5 


mini-spiral structure 















-8 





-8 


warped (inclined) disk? Need more information 



tical peak of (T») implies a broadening of the LOSVD in 
the bar region (Sect. 5.6). 

5.3. Embedded disks? 

In a significant fraction of the galaxies inner disks nested 
in the central regions of the large-scale bars are found 
from kinematical and morphological analysis (in NGC 470, 
NGC 2273, NGC 2681, NGC 5566, NGC 5905, and NGC 
6951) 6 . Below we summarize the signatures indicating this 
inner-disk structure: 



Here we do not mention NGC 3945 which has a massive in- 
ner disk embedded in the bar, discussed by Erwin et al. (2003). 
In his last paper based on photometric data Erwin (2004) also 
argued an inner disk origin for the isophotal twists in NGC 470 
and NGC 2273. 



— Photometry: an orientation of the inner isophotcs 
agrees with the PA of the line of nodes of an outer 
disk. 

— Kinematics: the central PAd yn of the stellar velocity 
field lies near the PA$. 

— Morphology: There are inner mini-spirals at distances 
corresponding to the nuclear disk (r < 1 kpc) in NGC 
2273, NGC 2681 NGC 5566, and NGC 6951. The ori- 
gin of the mini-spirals is not understood yet; various 
models are proposed (for review see Elmegreen et al., 
2002; Englmaier & Shlosman, 2000). Nevertheless, a 
"spiral" is some disk instability developed in the rela- 
tively thin and flat layer. The mini-spirals also appear 
in the images of NGC 4736 and NGC 7743, though we 
cannot determine inner disks in these galaxies. 

— Resonances: In the case of NGC 2681 and NGC 6951 
we have shown that an inner disk region is kinemat- 
ically decoupled because it lies inside the ILR of a 
global bar. 



18 



Moiseev et al.: Structure and kinematics of candidate double-barred galaxies. 



To confirm the inner disk hypothesis we tried to ex- 
tract their brightness distribution in the galactic plane. 
For this goal we subtracted a large-scale model of galaxy 
as a whole from the available images. A two-dimensional 
model of the surface brightness distribution is needed, be- 
cause barred galaxies are actually non-axisymmetric. Our 
model included a disk with central brightness /i^ and ex- 
ponential scale r<j, Sersic (1968) bulge and global bar. A 
two-dimensional surface distribution of the bulge bright- 
ness is characterized by three parameters: effective bright- 
ness n e , radius r e , and apparent ellipticity (see Eq. (1) 
in Moriondo et al., 1998). For a bar surface brightness 
distribution we adopted Ferrer's ellipsoidal distribution 
(Binney & Tremaine, 1987) projected onto the sky plane 
as: 

I(x,y)=I b (l-(x/af-(y/bfr (3) 

(the bar is aligned with the x-axis). Here a and b are 
the semi-axes of the bar, m is an integer parameter, lb (or 
fib when given in magnitudes) is the central brightness. 
Also the position angle of the bar PAb must be added for 
an arbitrary coordinate axes orientation. 

The model was calculated through the iterative tech- 
nique by means of the IDL-bascd software GIDRA 
(Moiseev, 1998). As a first step we average the original 
image data in the ellipses with PA and i and calcu- 
late a first approximation of the ma and Td in the disk- 
dominated region from the brightness profile obtained. 
The 2D model of the exponential disk is subtracted from 
the original image. The bar-dominated region of the resid- 
ual images then is interactively fitted for the determina- 
tion of the parameters of the bar. After removal of the 2D 
model of the bar we accept the ellipticity of the isophotes 
as a first approximation for the bulge ellipticity et and 
construct the averaged profile again. This profile is fit- 
ted by Sefsic's law. Then the cycle is repeated to refine 
the model. Usually 2-3 iterations are enough for a realis- 
tic model. The model-subtracted images and the accepted 
parameters of the models are shown in Fig. 21. We note 
that the profiles in the right plots of Fig. 21 arc calculated 
by averaging all component images, these are not a "clas- 
sical" one-dimensional fitting! Also we have deprojected 
the residual images to the galactic plane (Fig. 21, mid- 
dle column). The nuclear over-subtraction for NGC 470 
relates to the atmospheric seeing limitation. The depro- 
jected images reveal positive brightness residuals as true 
circular structures around the nucleus of NGC 470 and 
spiral structures in the inner round disks in NGC 2273 
and NGC 5566. The case of NGC 2681 is more complex, 
the residual disk is elliptical in the galactic plane, which 
can be connected with a turn of x 2 orbits inside the ILRs. 
Unfortunately there is no 2D kinematic information for 
r = 10 — 20", new observational data are needed for a 
more definite conclusion about this galaxy. 

The main conclusion is that the residual brightness dis- 
tribution shows a disk-like morphology (including ring-like 
one in NGC 470). Note that our decomposition has a pre- 
liminary character. More complex 2D models combining 



images in different bands (optical, NIR) and with differ- 
ent resolutions (HST plus ground-based) is required for a 
detailed inspection of the inner disks in barred galaxies. 
But such modelling is beyond the scope of our paper. 

5.4. Polar (inclined) disks 

The kinematic data provide evidence for inner polar disks 
in NGC 2681 and NGC 5850. Significant deviations (by 
more than 50°) of the PAd yn from the line of nodes occur- 
ring in the same direction as the PA of the inner isophotes 
are used as the main argument for this conclusion. Also a 
similar disk in NGC 3368 is detected on the basis of the 
HST morphology (the dust lane corresponds to the edge 
of the disk etc., see Sect. 4.6 and Sil'chenko et al., 2003) 
because the angular resolution of our 2D spectrographs 
fails at these spatial scales. What is the origin of these 
disks? We see two possibilities. 

Firstly, the disks could be formed by an accretion of ex- 
ternal matter with specific orientation of the angular mo- 
ment, as in "classical" large-scale stable polar rings (see, 
for example, Arnaboldi & Sparke, 1994). Indeed, NGC 
2681 demonstrates some spectral (widely distributed star- 
burst) and structural (disturbed dust lanes in the central 
disk) features, which have been ascribed to the merging of 
external gas (see Cappellari et al., 2001 and Sect. 4.4). It is 
possible that in NGC 3368 all ionized/neutral/molecular 
gas collected in the global warped disk has an external 
origin (Sect. 4.6 and Sil'chenko et al., 2003), which is is 
implied by the large-scale distribution of the neutral gas 
in the Leo I group (Schneider, 1989). Also the morphol- 
ogy of the neutral gas in NGC 5850 suggests a high-speed 
encounter with the nearby NGC 5846 (see Sect. 4.11 and 
Higdon et al., 1998). 

For the second point of view, we must note that the 
polar nuclear disks in NGC 3368 and NGC 5850 are or- 
thogonal to the major axes of the bars. In recent years, 
such polar mini-disks associated with a large-scale bar or a 
triaxial bulge have been detected in the circumnuclcar re- 
gions of several isolated galaxies, for example in NGC 2841 
(Afanasiev & Sil'chenko, 1999) or NGC 4548 (Sil'chenko, 
2002). It is possible that the gas in the centers of these 
galaxies has been moved into polar orbits due to the effect 
of the bar as is suggested in a qualitative model by Sofue 
& Wakamatsu (1994). Unfortunately, detailed numerical 
simulations of this process and a check of its stability are 
absent. 

It is possible that both mechanisms (internal as well 
as external ones) are at work in the described galaxies. 
For a detailed discussion on polar disks with sizes of a few 
hundred parsecs see Sil'chenko (2002) and Corsini et al. 
(2003a). In the last paper a list of 15 galaxies with nu- 
clear polar disks has been collected and their properties 
were investigated. New observations and numerical simu- 
lations are needed to understand completely these galactic 
structures. We hope that increasing amounts of observa- 
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Fig. 21. Surface brightness decomposition for four galaxies. Left column - the 2D model was removed from the original image. 
The solid line corresponds to the galactic major axis. Central column - deprojection of the residual brightness to the galactic 
plane. The gray scale (common for the left and central images) is in magnitudes. Circles mark the radius of the inner disks 
denoted in Table 8. Right column - surface brightness profiles. Solid lines: azimuthally averaged profiles of the images; dotted 
lines: the same for the components (bulge, bar, and exponential disk) of the model. The parameters of the components are also 
given here. 



tional information about inner polar disks will stimulate 
new theoretical work in this field. 

5.5. A counter-rotation 

We have found that the ionized gas in the central kilo- 
parsec of NGC 3945 rotates in the opposite direction with 
respect to the stellar disk. According to Kuijken et al. 
(1996), the gaseous disks in lenticular galaxies demon- 
strate a counter-rotation phenomenon in 24 ± 8% of all 
cases. Therefore, this is a frequent kinematic feature in 
early-type galaxies. This is probably attributable to a 
merger of an accreted gaseous cloud with the correspond- 
ing direction of angular momentum (Bertola et al., 1992). 
Against this background, counter-rotation of the gas in the 
center of NGC 3945, one of four SO galaxies in our sample 
(see Table 8) and one of three with gaseous disks (NGC 
2681, NGC 3945, and NGC 7743), comes as no surprise. 

5.6. Stellar velocity dispersion appearance 

Simple models of single-barred galaxies have shown (Miller 
& Smith, 1979; Vauterin & Dejonghe, 1997) that the cen- 
tral ellipsoidal peak in the distribution of must be 
aligned with the orientation of the outer bar. In our sam- 
ple the ellipsoidal peaks of the cr* distributions are found 
in six galaxies (NGC 470, NGC 2273, NGC 2681, NGC 
2950, NGC 3945, NGC 5905). As has been shown in our 
previous papers (Moiseev et at., 2002; Moiseev, 2002b), 
over our sample the position angle of the cr* peaks reveals 
a strong correlation with the PA of the major axis of the 
large-scale bar (see Fig. 4 in Moiseev, 2002b). But there 
is no any correlation between the symmetry axis of the 
velocity dispersion and the inner isophotal orientation. So 
in these objects the large-scale bars determine the dynam- 
ics of the stellar component even in the regions where the 
photometric "inner bars" are observed. 

In NGC 3786 a central drop of the stellar velocity 
dispersion is found (Sect. 4.7). Recently, Emsellem et al. 
(2001a) and Marquez et al. (2003) have observed similar 
drops in the 0% radial distributions in the central parts of 
some barred galaxies with AGN. These drops can proba- 
bly be attributed to the "cold" (from the dynamical point 
of view) stellar disk embedded into the bar. Indeed, the re- 
cent self-consistent N-body simulations by Wozniak et al. 
(2003) show that "young stars born in the nuclear regions 
from dynamically cold gas have a lower velocity disper- 
sion than the older stellar populations" . Such a drop was 
also previously observed in NGC 6951 (Perez et al., 2000; 
Marquez et al., 2003) up to distances of ±6 — 8" from 



the center (Fig. 5 in Marquez et al. , 2003) . Unfortunately, 
our cr* map for this galaxy has a smaller size (Fig. 12), so 
we could not detect this interesting feature of the stellar 
kinematics. 

5. 7. Molecular gas morphology 

Recently Petitpas & Wilson (2002, 2003) proposed a new 
method for resolving the origin of the secondary bar. They 
compared turns of NIR isophotes with molecular gas (CO) 
distributions. Namely, if NIR and CO images have simi- 
lar features, this may suggest that the isophote twists are 
related to structures hosted by disks and not to triaxial 
bulges. In our sample, we have found detailed CO-maps 
for NGC 470 (Jogee et al., 2001), NGC 2681 (Jogee et 
al., 2001), NGC 2273 (Petitpas & Wilson, 2002), NGC 
3368 (Sakamoto et al., 1999), NGC 4736 (Wong & Blitz, 
2000), and NGC 6951 (Kohno 1999). In five galaxies (ex- 
cept NGC 3368) the orientation of the CO-distribution 
coincides with the isophotes of our NIR (or optical) im- 
ages. We think that the physical reasons for this effect can 
be different: 

- In NGC 470, NGC 2273, and NGC 6951 the iso- 
contours of the inner CO distribution match the PA 
and e of the outer isophotes. We suggest that molec- 
ular gas is distributed in the nuclear disks discussed 
above (Sect. 5.3), as the radial scale of the CO concen- 
tration and the photometric disks are the same. Two 
gaseous clumps, symmetric with respect to the cen- 
ter found by Petitpas & Wilson (2002) in NGC 2273, 
can be connected with bright segments of the inner 
mini- spiral (Fig. 21), or may correspond to the stellar 
and gas orbit crowding near the ILR of the primary 
bar (see Petitpas & Wilson, 2002). Moreover, in the 
center of NGC 6951 the molecular gas distributed in 
the pseudo-ring agrees with the star-forming ring visi- 
ble in optical images (Sect. 4.13). But in contrast with 
the optical bands, the CO emission is dominated by 
a "twin peaks" morphology, which is also interpreted 
by Kohno (1999) as u xi/x 2 orbit-crowding regions". 
Therefore, the distribution of the molecular gas in the 
discussed galaxies can differ from purely disk-like, be- 
cause these disks nest in the non-axisymmetric poten- 
tial of a large-scale bar. Nevertheless, the CO emission 
is concentrated on the same scales as the photometric 
nuclear disks. 

- In NGC 2681 and NGC 3368 the central molecular 
gas is confined to the inner polar disks, hence these 
disks are CO-rich. But in NGC 3368 the NIR isophotes 
cannot be related to the polar disk, because this disk 
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is seen in the dust distribution only, not in that of the 
stellar component (see scheme in Fig. 18 in Sil'chenko 
et al., 2003). 

— In NGC 4736 there is a molecular bar present elon- 
gated along the stellar bar (see notes in Sect. 4.9). 

So the molecular gas morphology does not contradict 
the existence of the structures found in this paper (nuclear 
coplanar and polar disks). 

5.8. Where are the secondary bars? 

The circumnuclear regions (of a few kiloparsec in size) in 
all candidate double-barred galaxies of our sample (ex- 
cept NGC 5566) are dynamically decoupled, because de- 
viations of the gaseous/stellar dynamical axes relative to 
the lines of nodes or elliptical peaks of ct* are detected 
on our 2D maps. Various peculiar structures (inner po- 
lar and coplanar disks, mini-spirals) were revealed. Such a 
large fraction of peculiar structures is not surprising, be- 
cause different structural features can produce twists of 
the inner isophotes (see Introduction). However, we have 
not found any kinematic features connected with inner 
(secondary) bars. This cannot be explained by the limited 
angular resolution of our kinematic data, because we con- 
firm definitively the nuclear mini-bars (without large bars) 
in NGC 3368 and NGC 3786. We see two possibilities to 
explain this. 

Firstly, it is possible that the twist of the inner 
isophotes of galactic images has no connection with a dy- 
namically decoupled inner bar. In this case a "secondary 
bar" is an illusive structure, and different reasons may 
cause the isophotal distortion: inner disks or rings, nuclear 
spirals etc. Indeed, Shaw et al. (1993) had constructed a 
model where the inner twist of the stellar isophotes is trig- 
gered by a secular evolution of the gaseous component in 
a single-barred gas-rich galaxy. See also papers by Fricdli 
et al. (1996), Moiseev (2001b), and Petitpas & Wilson, 
(2002) where various mechanisms are considered for an 
isophotal twist explanation. 

Secondly, the kinematics of the inner bar in a double- 
barred galaxy differs essentially from single-bar dynam- 
ics. In this case, the deviations of PA^and the peaks 
of the velocity dispersion must be interpreted in a dif- 
ferent manner than in this paper. For example, in their 
last work Corsini et al. (2003b) have measured different 
pattern speeds for secondary and primary bar-like struc- 
tures in NGC 2950. However, PAd yn of the stellar velocity 
field coincides with the line of nodes (see Sect. 4.5). New 
self-consistent models of barred galaxies must be used to 
investigate this phenomenon. 

Of course, selection effects must be kept in mind be- 
cause our sample is small: we studied only 13 galaxies from 
Moiseev's (2001b) list, or only 6 objects from 50 "pho- 
tometrically confirmed" double-barred galaxies of Erwin 
(2004). However, our sample contains the most "classi- 
cal" double-barred objects with strong twist of the inner 
isophotes (NGC 2681, NGC 5850). 



Nevertheless, the existing samples of double-barred 
candidates must be revised. First of all, additional 2D 
kinematics data (panoramic spectroscopy, HI and CO in- 
terfcromctric radio observations) must be acquired for 
studying any dynamical decoupling of the circumnuclear 
regions. We hope that our paper is only the first small 
cautious step in this promising way. 

6. Conclusion 

The analysis of the morphology and kinematics of 13 
double-barred galaxy candidates observed by means of 
panoramic spectroscopy at the SAO RAS 6m telescope 
is presented in this paper. We tried to collect the maxi- 
mum possible information from observed 2D distributions 
of the optical and NIR surface brightness and kinematic 
properties. For this goal, panoramic spectroscopy observa- 
tions were carried out at the 6m telescope. Velocity fields 
and velocity dispersion fields of ionized gas and stars were 
constructed from these observations. The rotation curves 
and radial variations of the dynamical axis position angle 
were calculated. The analysis of the morphology of the 
galaxies was done by means of isophotc analysis of NIR 
and optical ground-based and HST archival images. 

We suggest that candidate double-barred galaxies are, 
in fact, galaxies with very different circumnuclear struc- 
ture. It is necessary to note that the majority of the ob- 
served morphological and kinematic features in our sam- 
ple galaxies may be explained without the secondary bar 
hypothesis. Three cases of inner polar disks, one counter- 
rotating gaseous disk and seven nuclear disks (with and 
without mini-spirals) nested in large-scale bars were found 
in this work. However, we think that the observational 
data presented here, particularly the 2D data, will be use- 
ful to proponents as well as to opponents of the secondary 
bar hypothesis. Doubtless the models of such galaxies 
must explain all their observed kinematic properties. 
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